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Abstract. We obtain a Central Limit Theorem for a general class of 
additive parameters (costs, observables) associated to three standard 
Euclidean algorithms, with optimal speed of convergence. We also pro- 
vide very precise asymptotic estimates and error terms for the mean 
and variance of such parameters. For costs that are lattice (including 
the number of steps), we go further and establish a Local Limit The- 
orem, with optimal speed of convergence. We view an algorithm as a 
dynamical system restricted to rational inputs, and combine tools im- 
ported from dynamics, such as transfer operators, with various other 
techniques: Dirichlet series. Perron's formula, quasi-powers theorems, 
and the saddle-point method. Such dynamical analyses had previously 
been used to perform the average-case analysis of algorithms. For the 
present (dynamical) analysis in distribution, we require estimates on 
transfer operators when a parameter varies along vertical lines in the 
complex plane. To prove them, we adapt techniques introduced recently 
by Dolgopyat in the context of continuous-time dynamics pTSJ. 



1. Introduction 

According to Knuth [201 p. 335], "we might call Euclid's method the grand- 
daddy of all algorithms, because it is the oldest nontrivial algorithm that 
has survived to the present day." Indeed, Euclid's algorithm is currently a 
basic building block of computer algebra systems and multi-precision arith- 
metic libraries, and, in many such applications, most of the time is spent 
in computing gcd's. However, the Euclidean algorithm has not yet been 
completely analyzed, and it is the purpose of this paper to provide such an 
analysis. 

We shall state informally our Central Limit Theorem (CLT) and Local Limit 
Theorem (LLT) later in this section, but first we discuss Euclidean algo- 
rithms and cost functions (equivalently: additive parameters) , recalling pre- 
viously known facts about their average-case and distributional analyses. 
Our results have been announced in 

Continued fraction expansions of real numbers. Every x g]0, 1] ad- 
mits a finite or infinite (Ci^)-continued fraction expansion of the form 

(LI) x= ^ — . 

mi H J 

m2 H J 

■■• + 

mn + ... 
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Ordinary continued fraction expansions can be viewed as trajectories of a 
one-dimensional dynamical system, the Gauss map T : [0, 1] [0, 1], 

T(x):=--|-J, for a; 7^0, TfO) = . 

X X 

(Here, [xj is the integer part of x.) For an irrational x, the trajectory 
T(x) = (x, r(x), r^(x), . . . , T"'(x), . . .) never meets and is encoded by 
the infinite sequence of digits {mi{x) , m2{x) , m^^x) , . . . ,mn{x), . . .), defined 
by mi{x) := m(T*~^(x)), with m{x) := [^J. It is usual to consider the 
truncated trajectory 7^(x) := (x,T(x), . . . ,T'^(x)), and let n tend to cxd. If 
X 7^ is rational, the trajectory T{x) reaches in a finite number of steps, 
and this number, -P(x), is called the depth of x. We set P{x) = oo for 
irrational x. 

Associate a nonnegative real value c(m) to each possible digit m > 1. We 
may also regard c as a function on the set 7i of inverse branches of T. We 
set for each x and each n 

min{n, P{x)) 

(1.2) Cn{x) := Yl c(m,(x)) . 

1=1 

We shall refer to c as a digit-cost and to Cn as the associated total cost (some- 
times, just cost), since a variety of costs, in the usual sense of computational 
complexity, are usually expressible in this way. 

The total cost (|1.'2() is a Birkhoff sum, i.e., a sum over iterates of the dynam- 
ics T. Hence, the functions c(m) play the role of "observables" in physics, 
and it is of interest to characterize their probabilistic behavior. Here, T 
admits a unique invariant density /i(x) = (l/log2)(l -|- x)~^, the Gauss 
density. When the integrability condition, 

(1.3, «c):=EcW-/ /.(.M.= rcW.))A(.M.<», 

Jh{[0,l]) Jo 

is satisfied, then the Birkhoff ergodic theorem implies that for Lebesgue 
almost every x in [0, 1] the average Cn{x)/n converges to /i(c) as n — > oo. 
We discuss next how transfer operators give refinements of this fact. 

The density transformer, also known as the Perron-Frobenius operator, 

(1.4) Hi[/](x) = ^|/i'(x)|-/o/j(x) 

was introduced early in the study of continued fractions (see e.g. Levy 
Khinchin [20], Kuzmin [HI], Wirsing [S], Babenko [Bl, and Mayer |31I). The 
density transformer is a special case of a transfer operator. Very general 
transfer operators were introduced by Ruelle, in connection with his ther- 
modynamic formalism (see e.g. IH|). We shall see that "weighted" transfer 
operators appear naturally in the probabilistic analysis of dynamics. 

Fix a reference probability measure on [0,1], absolutely continuous with a 
smooth density /, and denote by E[-] the corresponding expectation. To 
establish probabilistic results on total costs along truncated trajectories as 
n ^ CXD, it is standard to use the sequence of moment generating functions, 
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i.e., the sequence of expectations E[exp{wCn)], for w complex. In probabilis- 
tic i.i.d. situations, E[exp(t(;C„)] is the n-th power of some expectation. In 
our setting, a quasi-powers approximation (with remainder term) can be ob- 
tained for E[exp(tt;C„)] after expressing it in terms of the nth iterate H"^ 
of a transfer operator. The transfer operator in question is the following 
perturbation of the density transformer: 

(1.5) Hi,^[/] = J^expM/i)] • \h'\ -foh. 

heh 

The density transformer Hi = Hi^o (acting on functions) has a dominant 
eigenvalue A = 1, and a spectral gap: the rest of the spectrum lies in a 
disk of radius < 1. For a non constant digit-cost satisfying a moderate 
growth condition (|2.5j) . elementary perturbation theory |28j implies that Hi^t„ 
inherits the spectral gap when w is near 0. This gives the above-mentioned 
quasi-powers expansion and, together with convexity of the logarithm of the 
dominant eigenvalue, leads to a proof that the asymptotic distribution of 
the total cost is Gaussian, with (optimal) speed of convergence 0{l/^/n). 
This Central Limit Theorem [stated below more precisely as Theorem 1] 
is quite well-known. See for instance ^Hl or [2] for interval maps, and |^ 
for a more abstract framework and references to the pioneering paper of 
Nagaev. It is convenient to base the proof below on a compact and versatile 
statement, the "Quasi-Powers Theorem" of Hwang (^I^HIUZI [Theorem 
below] , which encapsulates the consequences of the Levy continuity theorem 
and the Berry-Esseen inequality. 

Continued fractions of rational numbers: Average-case analysis 
of Euclidean algorithms. There are variants of the standard continued 
fraction algorithm induced by variations of the standard division procedure. 
See |49j for many examples and a classification into the "Fast" and the 
"Slow" Class. In this paper, we study three algorithms of the Fast Class, 
the standard, centered, and odd algorithms, specified at the beginning of 
Section 2. An execution of a Euclidean algorithm on the input {u, v) formed 
with two integers u, v such that u/v = x gives rise to a rational trajectory 
T(x) which ends at zero in P{x) steps, and the total cost of the trajectory 
is then 

P{x} 

(1.6) C{x) := c(m,(x)) . 

i=l 

The reference parameter is no longer the truncation degree n, but the 
size N := msix{u,v) of the input {u,v). The reference probability mea- 
sure Pat is now the uniform discrete measure on the (finite) set of inputs of 
size < N. 

As it is often the case, the discrete problem is more difficult than its con- 
tinuous counterpart. However, the average-case complexity of Euclidean al- 
gorithms is already well-understood, as we explain next. For the number of 
steps P{u, v), which corresponds to the trivial digit-cost c = 1, the standard 
Euclidean algorithm was first analyzed in the average-case around 1969, in- 
dependently by Heilbronn ,22, and Dixon [15| . The centered algorithm was 
studied by Rieger [5^ . 
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Consider now a general digit-cost c of moderate growth and the associated 
total cost C of the rational trajectories. The expectation EAr[C] is described 
by the partial sums of the coefficients of a generating function S{s) (a com- 
mon tool in the average-case study of algorithms jl9l l2Uj ^ where the pa- 
rameter s "marks" the size N of inputs. As it is usual in number theory, 
the generating functions S{s) are Dirichlet series. Recently, Vallee [4^] has 
related S{2s) to the quasi-inverse (/ — H^)^^ of another perturbation of 
the density transformer, together with its weighted version : 

(1.7) H,[/] = \h'\^ -foh, := c{h) • \h'\^ -foh. 

hen hen 

Then, spectral information on may be used to show that (/ — H^)^^ is 
analytic in the half-plane {?R.s > 1}, and analytic on 5Rs = 1 except for a 
simple pole at s = 1. Under these conditions, one can extract asymptotically 
the coefficients of S{s) by means of Delange's Tauberian theorems |14[ I44j. 
(Just like in the usual prime number theorem, or in weighted dynamical 
prime number theorems, see e.g. Sections 6-7 of |37j.) For costs of mod- 
erate growth and Euclidean algorithms in the Fast Class, this dynamical 
approach gives that the mean value E7v[C] of the total cost of the ratio- 
nal trajectory satisfies IE7v[C] ~ /1(c) • filogN. Here, /i(c) is the asymptotic 
mean value (|1.3j) of truncated real trajectories, and /i equals 2/|A'(l)|, where 
A(s) is the dominating eigenvalue of Hg. 

Euclidean algorithms: Distributional analysis and dynamical 
methods. Main results. We have seen that, with respect to any cost 
of moderate growth, rational trajectories behave in average similarly to the 
way truncated real trajectories behave almost everywhere. It is then natu- 
ral to ask whether this analogy extends to distributions: Is it true that the 
distribution of the total cost C{x) on rational trajectories with an input x 
whose numerator and denominator are less than is asymptotically Gauss- 
ian (when N tends to oo)? How to compare the distribution of some cost on 
truncated real trajectories and on rational trajectories? This paper provides 
a precise answer to all these questions for three different algorithms that all 
belong to the Fast Class. 

Concerning the standard Euclidean algorithm and the number of steps (i.e., 
the constant cost c = 1), Hensley ^24^ has obtained a Central Limit Theorem, 
and a Local Limit Theorem with speed of convergence 0((log A)^^/^^). In 
the present work, we apply dynamical methods for the first time to the dis- 
tributional analysis of discrete algorithms; in this way, we improve Hensley's 
result while extending it to a large class of cost functionals and to several 
algorithms. 

Our strategy consists in describing the moment generating function 
Ejv[exp(t(;C)] as a quasi-power. It turns out that KN[exp{wC)] is related to 
the partial sums of the coefficients in a bivariate series S{s,w). This series 
is of Dirichlet type with respect to the variable s, while the extra param- 
eter w "marks" the cost c, and we require uniform estimates with respect 
to w. Tauberian theorems are now insufficient, since they do not provide 
remainder terms: We need a more precise "extractor" of coefficients, and 
the Perron formula is well-suited to this purpose. 
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We establish below (|2.17() a simple relation between S{2s, w) and the quasi- 
inverse of a two-variable transfer operator 'H.s,w defined by 

Hs,«,[/] = ^ expM/i)] • |/i'r / o . 
hen 

Note that this operator is a simultaneous extension of the three operators 
Hi, Hi^it,, Hs defined in (|1.4|) . (|1.7|) . and while ni^^ is just its derivative 
with respect to w at w = 0. In order to apply Perron's formula with a 
convenient integration contour, it is thus useful to know that, in a half 
plane "perforated" at s = 1, of the form {5Rs>l — e,|s — 1|> e/2} for small 
e > 0, a certain norm of the quasi-inverse satisfies 

(1.8) ||(/-H,,^)-i||, < max(l,|9s|«), 

with ^ < 1, uniformly in w. Then, spectral properties of H^^^ inherited 
from Hi give the desired quasi-power expansion for E,]\f[exp{wC)], for w 
close to 0, and some Pat close to Pat. Note that Hensley in ]^ has used a 
transfer operator Hg^O; but in an appreciably different way: Hensley obtains 
distributional results on rational trajectories upon aproximating discrete 
measures on rationals by continuous measures. In particular, his approach 
avoids parameters s of large imaginary parts. 

Adapting powerful methods due to Dolgopyat |12J, we show [Theorem 2 
and Lemma 6 below] that the quasi-inverse satisfies the estimates (|1.8|) for 
large \Qs\. Dolgopyat was interested in the decay of correlations for hyper- 
bolic flows satisfying some uniform nonintegrability condition ( UNI) . Later 
on, Pollicott and Sharp used Dolgopyat's bounds together with Perron's 
formula to find error terms in asymptotic estimates for geodesic flows on 
surfaces of variable negative curvature; see e.g. where only univariate 
Dirichlet series with positive cofficients appear. To the best of our knowl- 
edge, the present paper is the first instance where these powerful tools are 
extended to dynamical systems with infinitely many branches and applied 
to distributional analyses in discrete combinatorics. 

Let us now state informally our two main results about the three algorithms: 
Theorem 3. Consider a nonzero digit-cost c of moderate growth. We show 
the following Central Limit Theorem: the asymptotic distribution of the 
total cost C{u/v) of an execution of the algorithm on the rational input 
u/v, uniformly randomly drawn from {{u^v) ,u < v < A^}, is asymptotically 
Gaussian, with best possible speed of convergence, of order 0{1/ \J\ogN). 
We give expansions for the expectation and variance, which are asymptoti- 
cally proportional to log A. The constants /x(c) and 5^(c) in the main terms 
of the expectation and the variance are expressed in function of the partial 
derivatives at {s,w) = (1,0) of the dominant eigenvalue of tis^w, and alter- 
natively in terms of /i = /i(l), (5^ = (5^(1), and the constants iu(c), 5{c) from 
Theorem 1. In particular /x(c) admits a closed form. 

Theorem 4. A digit-cost is lattice if it is non zero and there exists L G 
so that c/L is integer- valued. The largest such L is then called the span of c. 
For instance, any nonzero constant c is lattice. For lattice costs of moderate 
growth, we obtain a Local Limit Theorem with optimal speed of convergence 
0(l/i/log A). This time, we use estimates for E7v[exp(irC)], where r varies 
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in a compact set of the real line. They lead, with the saddle-point method, 
to a very natural and concise proof. 

Three special instances of our results for lattice costs are of major interest. 

(i) Digit-cost c = 1. For each of our three algorithms the number of steps is 
asymptotically Gaussian, with mean fi log (where fi admits a closed form) 
and variance 6'^ log N, in the sense of the CLT and the LLT with speed 

o(i/yb^). 

(a) Digit-cost c = Cm = the characteristic function of a digit m. The 
total number of occurrences of a fixed digit m in a rational trajectory for 
our three algorithms is asymptotically Gaussian (CLT, LLT), with mean 
yu • /I(cm)logA^. The constant Ji{cm) is explicit for the three algorithms. In 
the case of the standard Euclidean algorithm we recover (see [3^]) 

/^(Cm) = log (l + —^^—rr) . 

log z \ m(m + 2) / 

{Hi) Digit-cost c = the binary length I of the digit. The binary encoding 
of a rational trajectory is asymptotically Gaussian (CLT, LLT), with mean- 
value /i • '{l{t)\ogN . The constant /I(^), explicit for the three algorithms, 
is a variant of the Khinchine constant |29j . For the standard Euclidean 
algorithm, it is equal to (see [l9] ) 

oo 
fe=0 

Plan. In Section 2, after describing our three Euclidean algorithms and 
their underlying dynamical systems, we state and prove Theorem 1 (using 
Theorem 0) and state our main results. Theorems 2 (our version of Dolgo- 
pyat's bounds), 3 (CLT), and 4 (LLT). The first part of Section 3 is devoted 
to proving Theorem 2 and checking that its assumptions hold for our three 
algorithms. In the last two subsections of Section 3, we obtain quasiperi- 
odicity results. With Theorem 2, they entail expected properties for the 
Dirichlet series. In Sections 4 and 5, we present the proofs of Theorems 3 
and 4, respectively. 

2. Dynamical methods and statement of results. 

After a description of the three Euclidean algorithms to be studied and their 
associated dynamical systems (§2.1), we introduce the weighted transfer 
operator in §2.2. Section 2.3 explains the fundamental role this operator 
plays in the distributional analysis of truncated real trajectories. Theorem 1 
is stated there, and proved using Hwang's Quasi- Power result (Theorem 0). 
Next, we turn to rational inputs and Euclidean algorithms. We introduce in 
§2.4 Dirichlet series of moment generating functions, which we relate to the 
quasi-inverse of the weighted transfer operator. We briefly explain in §2.5 
how to apply the Perron formula to extract coefficients of Dirichlet series, 
stating also estimates a la Dolgopyat (Theorem 2) useful for the applicability 
of this formula. Finally, in §2.6 we state our CLT (Theorem 3) and our LLT 
for lattice costs (Theorem 4). 
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2.1. Euclidean algorithms and their associated interval maps. 

Three Euclidean algorithms are to be analyzed; each of them is related 
to a Euclidean division. Let v > u > 1 he integers. The classical division, 
corresponding to the standard Euclidean algorithm Q, v = mu + r produces 
an integer m > 1 and an integer remainder r such that < r < u. The 
centered division (centered algorithm /C) requires v > 2u and takes the form 
V = mu+s, with integer s G [—u/2, +u/2[. Letting s = er, with e = ±1 (and 
e = +1 if s = 0), it produces an integer remainder r such that < r < u/2, 
and an integer m>2. The odd division (odd algorithm O) produces an odd 
quotient: it is of the form v = mu + s with m odd and integer s G [—u, +u[. 
Letting s = er, with e = ±1 (and e = +1 if s = 0), it produces an integer 
remainder r with < r < u, and an odd integer m > 1. In the three cases, 
the divisions are defined by pairs q = (m, e), which are called the digits. 
(See Figure 1.) 

Instead of the integer pair (u, v), we may consider the rational u/v, since the 
pair (du, dv) produces the same sequence of digits as {u, v), up to multiplying 
all remainders r by d. Then, the division expressing the pair [u, v) as a 
function of (r, u) is replaced by a linear fractional transformation (LFT) h 
that expresses the rational u/v as a function of r/u. For each algorithm, 
the rational u/v belongs to the interval T' defined in Fig 1. 

To summarize, on the input {u,v), each algorithm performs a sequence of 
admissible Euclidean divisions, of the form v = mu + er with r/u & X'U {0}, 
and (m, e) G Di. On the input u/v, it performs a sequence of LFT's from a 
generic set 7i (depending on the algorithm) whose elements h^^j are indexed 
by the admissible pairs q = (m, e) of Pi and are of the form h[jn^^]{x) = 
l/(m + €x). The LFTs appearing in the final step belong to a subset T dTi 
and are indexed by the admissible pairs q = {m, e) of T>i H T>2- (Fig. 1). 

Each algorithm applied to a rational u/v builds a specific continued fraction 
(2.1) 



mi + 



m2 + 



€2 



mp 



of depth P, decomposing also u/v as 

(2.2) {u/v) = /ti o /i2 o . . . o /ip(0) = /i(0) , 

with hi ^ 'H, \ < i < P — \, and hp e J^. We are interested in various costs 
related to an execution of the algorithm. The most basic one is the number 
of steps P{u,v). In general, given a digit-cost function c on the set H, we 
consider additive cost functions C of the form 

P{u,v) 

(2.3) Ciu,v) := J2 

i=l 

We next see how to associate to each algorithm a dynamical system of the 
interval, T : T ^ T. The interval X is defined in Figure 1. The map 
T extends the map defined on rationals by T{u/v) = r/u, where r is the 
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Algorithm 



Q (standard) 



K (centered) 



O (odd) 



Intervals 



J =[0,1], T'=]0,1[ 



J =[0,1/2], J'=]0,l/2] 



J =[0,1], J'=]0,1] 



Generic conditions 
Vi on pairs (m, e) 



m > 1, e = +1 



m> 2, e = ±1 

if m = 2 then e = +1 



m > 1 odd, e = ±1 
if m = 1 then e = +1 



Final cond. 



m > 2 



e = +l 



e = +l 



ililllH ! i 'i 



Graph of 

n^) = \-.-M^)\ 



Function A{y) 



Integer part of y 



Nearest integer to y 

i.e. m s.t. y-me [-1/2, +l/2[ 



Nearest odd integer to y 

i.e. m odd s.t. y — m £ [—1, +1[ 



Gontraction ratio 



p= 1/(72 + 1)2 



p=l/</.2 



~T 1 — 

log2 1 + x 



1 , i 



T" 



T" 



T" 



Invariant density fi 



< + X (jr — X 



+ 



X (jr — X 



Entropy 



6 log 2 



61og(^ 



Qloggji 



Figure 1. The three Eudidean systems (^ = (1 + \/b)/2). 



remainder of the Eudidean division on (u,v). We get 



T{x) :- 











X \ 





, x^O, T(0)=0, 



A depends on the algorithm and is defined in Figure 1. It is easy to see that 
the maps associated to the three algorithms belong to the following class: 

Definition. [Piecewise complete maps of the interval] A map T : I ^ I is 
piecewise complete if there exist a (Gnite or countable) set Q, whose elements 
are called digits, and a partition {Iq}qeQ (modulo a countable set) of the 
interval I into open subintervals Ig such that the restriction of T to Ig 
extends to a a bijective mapping of class from the closure of Tg to X. 

We may consider the set H. = {h^g]} of branches of the inverse function T^^ 
of a general piecewise complete map T, naturally indexed by the set Q. The 
set of the inverse branches of the iterate is 7i^; its elements are of the 
form h[q^j o h^^^j o • • • o h^^^j where k is called the depth of the branch. Setting 
= {Id}, the set H* := Dk>o>i^ is the semi-group generated by Ji. Each 
interval h{T) for h of depth k is called a fundamental interval of depth k. 
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2.2. Transfer operators of interval maps. The maps (X, T) associated 
to the three Euchdean dynamical systems belong in fact to a subclass of 
piecewise complete mappings : 

Definition [Good class] . A piecewise complete interval map {I, T) belongs 
to the good class if: 

(i) T is piecewise uniformly expanding, i.e., there are C and p < 1 so that 
\h'{x)\ < CfP for every inverse branch h of T", all n and all x £ I. The 
infimum of such p is called the contraction ratio, and satisfies 

(2.4) p = limsup(max{|/i'(x)|;/iG J})^/" . 

n— >oo 

[ii) There is K > called the distortion constant, so that every inverse 
branch hofT satisfies \h"{x)\ < K\h'[x)\ for all x ^Z. 

{in) There is ctq < 1 such that 'YlheH \ < °° -^^^ ^^^^ o" > ctq. 

(Note that maps in the good class are automatically topologically mixing.) 
The distortion condition in {ii) follows from the assumption when there 
are finitely many branches, and is just Renyi's condition otherwise. 

To check that maps associated to our algorithms are in the good class, for 
(To = 1/2, use = 0(m^^). (See also Figure 1, and [12] for proofs.) 

If T is endowed with an initial probability density with respect to Lebesgue 
measure, T acts on it and transforms it into a new density gi. The operator 
H such that gi = ^[go] is called the density transformer, or the Perron- 
Probenius operator (acting now on functions, soon we shall restrict its 
domain) . An application of the change of variable formula gives 

H[/](x) ■.= Y,\h'{x)\ foh{x). 

hen 

It is useful to deal with a more general operator, the transfer operator Hg 
which depends on a complex parameter s: 

H,[/](x):=^|/i'(x)r ■foh{x). 

hen 

(Note that Hi = H.) If cr := 5Rs > do, then acts boundedly on the 
Banach space C^(T) of functions on X endowed with the norm 

= ll/llo + ll/lli, with ll/llo := sup I/I, ll/lli := sup |/'|. 

We consider nonnegative cost functions satisfying the following condition: 

Definition [Condition MQ -Moderate growth-]. Let Ti. be the inverse 
branches of a map in the good class. A digit-cost c : TC ^ M'^ is of moderate 
growth if it is non identically zero and if the series 

(2.5) ^exp[u;c(/i)] • |/i'(x)|" 

hen 

converges when (Ks, "Siw) belongs to a real neighborhood of (1,0), or, equiv- 
alently to Sq x Wq with Sq =]S'o, +oo], for ao <ao < 1, and Wq =] — oo, i^qI 
for vq > 0. 
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Definition [Class GAiG] A triple (X, T, c) formed with an interval map 
T : I ^ Z of the good class, and a digit-cost c of moderate growth will be 
said to be of QMQ-type. 

We extend the digit-cost to a cost function, also denoted c, on Ti* by 

k 

c{hi o /i2 o . . . o /ifc) := ^ c{hi) . 

1=1 

We can now define the weighted composition operator which depends on 
two (complex) parameters s and w, 

(2.6) H,,^[/](x) := expM/i)] • \h'{x)\' ■ f o h{x) . 

The additive property of costs and the multiplicative property of the deriva- 
tives entail 

Hr,^[/](x) := ^MMh)] ■ \h'{xr ■ f o h{x), 

and the quasi-inverse (/ — Hs^^„)~^ can be written (formally) as 

(2.7) (/ - H,,^)-i[/](x) := Y exp[u;c(/i)] • \h'ix)\' ■ f o h{x) . 

hen* 

We recall next some well-known spectral properties of the transfer operator 
tls,w Endow the Banach space C^(I) with the norm || • It is known 
that for (JRs, Rw) G Sq x Wq the operator H^^^ is bounded but not compact 
acting on C^iT) [see e.g. ^JIH]]; however, it is quasi-compact. We recall the 
definition of quasi- compactness for a bounded operator L on a Banach space: 
Denote by Sp L the spectrum of L, by -R(L) its spectral radius, and by i?e(L) 
its essential spectral radius, i.e., the smallest r > such that any A G Sp(L) 
with modulus |A| > r is an isolated eigenvalue of finite multiplicity. An 
operator L is quasi-compact if -Re(L) < -R(L) holds. 

We denote the partial derivatives of first and second order of a function 
F{s,w) at {a,h) by F;(a,6), F',{a,h), F^,{a,b), F'^,{a,h), F::„{a,b). 

Proposition [Classical spectral properties of transfer operators]. Let 
ils,w be the transfer operator i)Z6'|) associated to a GM.Q triple (T, T, c) with 
contraction constant p. Denote by R{s,w) its spectral radius and Re{s,w) 
its essential spectral radius. Let Sq, Wq be the real sets from i)2.5|) . When 
w = 0, we omit the second index in the operator and its associated objects. 

(1) [Quasi-compactness.] Let p < p < 1. If {a = ^s,u = ^w) G Sq x 
Wq, then iis,w d-cts boundedly on C^{T). Then R{s,w) < R{Rs,Rw) and 
Re{s,w) < p ■ R(Jfis,?fiw), in particular iis,w is (uniformly) quasi-compact 
for real {s,w). 

(2) [Unique dominant eigenvalue.] For real {a,^) G Sq x Wq, Ho-^jy has 
a unique eigenvalue X{a,i') of maximal modulus, which is real and simple, 
the dominant eigenvalue. The associated eigenfunction f^^y is strictly pos- 
itive, and the associated eigenvector p^j^y of the adjoint operator H* is a 
positive Radon measure. With the normalization conditions, /icr,;/[l] = 1 
and Pa,u[fa,u] = 1, the measure pa,u '■= fa,uV-a,u is a probability measure. In 
particular, pi is Lebesgue measure, with A(l) = 1. 
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(3) [Spectral gap.] For real parameters {(t,v) G Sq x Wq, there is a 
spectral gap, i.e., the subdominant spectral radius r^^y > Re{c^, defined by 
r„^y := sup{|A|;A e Sp(H^^,,), A / A(fT, i/)}, satisfies r^^^ < X{a,i'). 

(4) [Analyticity in compact sets.] The operator Hs^w depends analyt- 
ically on {s,w) for {^s,^w) G Sq x Wq. Thus, A((T, z^)^-^, f^l, and 
depend analytically on (a, i^) G Sq x Wq, and are uniformly bounded in any 
compact subset. 

(5) [Analyticity in a neighborhood of (1,0).] If {s,w) is complex near 
(1,0) tiien X{s,w)^^, f^^, and /^^^ are well-defined and analytic; moreover, 
for any 6, with ri < 9 < 1, one has ri^w/\^i'^,w)\ < 9. 

(6) [Derivatives of the pressure.] For (cr, v) £ T^qX Wq, define the pressure 
function A{a, v) = log A (cr, u). 

(6. a) A'(l) is the opposite of the Kolmogorov entropy of the dynamical 
system (T, ^ui). Also, A'^(1,0) is the /ii -average of the cost: 



A'(l) = - / log|r'(x)|/i(x)dx<0, A;(1,0) = Vc(M / fiix)dx. 



(G.b) If c is not constant, the second derivative A^2(l, 0) is strictly positive. 

(7) [Function w (^(w).] There is a complex neighborhood W of and 
a unique function a : W ^ C such that X{a{w),w) = 1, this function is 
analytic, and o"(0) = 1. 

Proof. We refer to [IJ EH] except for claim (7) (which follows from 
A'(l) 7^ and the implicit function theorem) and for claim (6): 
(6. a) Taking the derivatives at (1,0) oiiis^w[fs,w] = ^{s,w)fs^w (with respect 
to s or w), integrating on I with respect to /ii^o (equal to the Lebesgue 
measure), and using that preserves Jli, gives the expressions as integrals. 
To finish, apply Rohlin's formula. 

(6.6) Convexity of the pressure is an old theme, see, e.g., [IJ and also 
PTTj . Chapter 4, Prop. 10 in gSI, Prop. 3.8 in Proposition 6.1 in ^. 
We adapt here the last work. It is stated in the context of functions with 
bounded variation. Due to our strong Markov assumption, we may work in 
C^{1). Since /i is a strictly positive function and Hi[c] belongs to C^, 
we may transfer Broise's proof to our context: it shows that A|^2(l,0) is 
zero if and only if there exists a constant K and a function u £ C^(T) for 
which, for all h £ 7i, the equality c{h) = u — u o h + K holds. Using the 
fixed point of each branch h proves that c is constant. □ 

We have already remarked that the three Euclidean dynamical systems be- 
long to the good class. A sufficient condition for the cost c to be of moderate 
growth is c(m, e) := c(/i[m,e]) = O(logm). Note that, for an inverse branch 
of depth k, of the form /ijg^] o • • • o /ijg^j, the interval h{I) gathers all real x 
for which the k first digits of the CF expansion are {qi,q2, . . . ,qk)- Further- 
more, each inverse branch of any depth is a linear fractional transformation 
h{x) = {ax -\- h)/{cx + d), with a, b, c, d coprime integers, with determinant 
ad — be = ±1, and denominator D[h] related to \h'\ through: 




D[h]{x) := \cx + d\ = \deth\ 



/^\h'{x)r'/^ = \h'{x)\-'/\ 
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Therefore, the transfer operator can be alternatively defined by 

(2.8) H,,^[/](x) = expK(M] ^ ^^ f{h{x)) . 

The above reformulation will be useful in §2.4. 

We associate to each of the three algorithms and its final set JT c a 
transfer operator 

(2.9) ^ 

Fs,yo[f]ix) := exp[wc{h)] f{h{x)) = Hs,^[/ • lu^^:^h{i)]ix). 

her V \\ I 

It is easy to see that Fs,«, acts boundedly on C^(X) for (5Rs, 5Rt(;) E Sq x Wq, 
and to generalize the relevant statements of Proposition to this operator. 

2.3. Transfer operators and real trajectories. We consider one of our 
three algorithms, or more generally any triple of QMQ type, with a non- 
constant cost c. The interval I is endowed with a probability measure with 
smooth density / and the set of endpoints of fundamental intervals (rational 
points for Euclidean algorithms) can be neglected. We are interested in the 
asymptotic behavior of the distribution of C„(x) := Yl'i=i'^i^i) when the 
truncation degree n tends to oo. As already mentioned, a very convenient 
tool is the Levy moment generating function of the cost, E[exp(wC„)] = 
YlheTf^ exp[zi;c(/i)] • j^^j^ f{y) dy. The change of variables y = h{u) gives 
(2.10) 

E[exp{wCn)] = [ Y expK(/i)] • \h'{u)\ ■ f o h{u) du = j Yll^^[f]{u) du . 
■^^ hen" 

The above relation is fundamental for analysing costs on truncated real 
trajectories, as we explain next. 

By Proposition 0, for any 9 with ri < 9 < 1, there is a small complex 
neighborhood W of 0, so that, for w G W, the operator Hi^„, splits as 
Hi^^ = X{l,'w)Pi^w + ^i,w, where Pi,^ is the projector for the eigenvalue 
X{l,w) and R(Ni^w) < 9X{l,w). Therefore, 

HlJfKu) = A(l,i«)"Pi,^[/](n) + NlJf]iu) , Vn > 1 , 

which entails 

E[exp(u;C„)] = (^A(l, w;)" Pi,^[/](n) d^x^ (l + 0(r)) , 

with a uniform O-term for w € W. In other words the moment generating 
function behaves as a "quasi-power," and we may apply the following result: 

Theorem 0. [Hwang's Quasi- Power Theorem] |251l26| l27j Assume that the 
moment generating functions for a sequence of functions Cm on probabihty 
spaces {il,N,¥]\f) are analytic in a complex neighborhood W of zero, where 

(2.11) ENiexpiwCN)] = expiPNUiw) + Viw)] (l + 0(k^^)) , 

with f3j\f, Kj\f ^ oo as N ^ oo, and U{w), V{w) analytic on W. Assume 
C/"(0) 7^ 0. Then, the mean and the variance satisfy 

En[Cn] = PnU'{0) + y (0) + 0(k^1) , 
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Viv[Cjv] = PnU"{0) + V"{0) + 0(«;^^) . 

Furthermore, the distribution ofCNonQ,^ is asymptotically Gaussian, with 
speed of convergence 0{k^ + 

For each fixed k > 3, there is a polynomial of degree at most k, with 
coefficients depending on the derivatives of order at most k at Q of U and 
V, so that the moment of order k satisfies 

(2.12) En[C^] = PkiPN) + O 

with a O—term uniform in k. 

Proof. [Sketch.] This statement encapsulates a classical calculation analo- 
gous to the proof of the central limit theorem by characteristic functions.The 

speed of convergence then results from the Berry-Esseen inequalities, the 
moment estimates are consequences of the derivability of analytic func- 
tions. □ 

For our application, we set 0,n = i^:fdx) for all n, Cn = Cn, Pn = n, 
Kn = 9"'^. The function U is the pressure function w i— A{l,w), and 
V{w) = log (^Jj-Pi^yj[f]{u) du) . Since c is not constant, the function A{l,w) 
is absolutely convex at w = 0, (see Proposition 0) and U"{0) / 0. Thus, 
using the formula for A!^(1,0) in Proposition 0, Theorem entails the fol- 
lowing Gaussian asymptotic distribution result, which applies in particular 
to our three Euclidean algorithms. 




Theorem 1. For a triple {I,T,c) of QAiQ type with non-constant c and 
any probability F on I with a density, there are /i(c) > and 6{c) > 
so that for any n, and any F G M 



Cn{x) - fl{c)n 
5{c)^/n 



< Y 



Furthermore, (recalling that vi is the subdominant spectral radius of the 
density transformer Hj, for any 9 which satisfies r\ < 9 <1, one has : 



E [Cn] = Kc) ■ n + r7(c) + 0{9^) , V [C„] = ? (c) • n + <5i (c) + O(0") , 



with ji{c) = a;(i,o) = a;(i,o), P{c) = a:.(i,o) = a" (i,o) - a;^(i,o). 



Finally, fi{c) involves the invariant density f\ , 

= y]c(g) / h{x)dx. 



Note that /i is explicitly given in Figure 1 for Q, O, and /C, so that /2(c) is 
computable in these cases. 
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2.4. Dirichlet generating functions and transfer operators. We re- 
strict now our study, for each of our three algorithms, to (nonzero) rational 
inputs X = {u/v) S X' and a cost c of moderate growth. The intervals X' are 
defined in Figure 1. We consider the sets 

O := {(u, e , - G J'} , VI := {(u, u) E ; gcd (u, = 1} , 

V 

and we endow the sets 

VLn ■■= {(m, v) £Q.;v < N} , Q.N := {(u, v) £^;v < N}} , 

with uniform probabilities Pat and P^v, respectively. For the moment we 
only consider 17 at. 

To study the distribution of the total cost C{u,v) 1)2. 3() associated to some 
digit-cost c and restricted to Oat we use its moment generating function: 

(2.13) EA.[exp(«;C)]:=|^, 

where <^>^(A^) = (^c,w[N) is the cumulative value of exp[it;C] on 17 at, 

(2.14) ^u,{N):= Yl eMwC{u,v)], ^N) = \Qn\ ■ 

(M,ti)enjv 

Extending the principles defined in |^ UHl OHj , we replace the sequence of 
moment generating functions by a single bivariate Dirichlet series, hence- 
forth called the Dirichlet-moment generating function: 

(2.15) S{s,w):= ^^M^C{u,v)]=Y^^ 

{u,v)£fl n>l 

where Cniw) := J2{u,v)en„ ,v=n(^MwC{u,v)]. 

Since the partial sum of the coefficients of the series S{s,w) satisfies 

(2.16) c„(?i;) = $^(7V), 

n<N 

to analyze the moment generating function EAr[exp(wC)] of the cost C on 
riAT, it suffices to estimate the functions $^(A^) (asymptotically in — > cxd, 
and uniformly in it; in a complex neighborhood of 0). 

As we previously did for truncated real trajectories, we aim to relate the 
moment generating function of costs on rational trajectories to the weighted 
transfer operator. An execution of the Euclidean algorithm on a input 
(u, v) S performing P{u, v) steps uniquely decomposes the rational 

u 

— = /ii o /i2 o . . . o hp(0) = h(0), 

V 

with hi £ 7i,l < i < P — 1 and hp £ J-. Thus, each Euclidean algorithm 
defines a bijection between the sets and TL* x F. In view of 1)2. 8112. 9() and 
1)2. 15() . the relations 

P{u,v) 

v = D[h]{{)), C7(u,i;) = c(/i) = ^ c(/ii) 

1=1 
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provide the desired expression for the Dirichlet moment generating function 
S{s,w) in terms of the transfer operators Hs^w and Fg^w- 

(2.17) 5(2s,ti;) =F,,^o(/-H,,^)-^[l](0). 

Returning to il, we remark that each element {u',v') of Q can be written 
in a unique way as {du, dv) with d = gcd(M, v) and (n, v) G VL, and, as 
already observed, an execution of the algorithm on (n', v') leads to the same 
decomposition p.2|) . the two costs C{u',v') and C{u,v) being equal. We 
may define a Dirichlet moment generating function S{s,'w), and we get 
(2.18) 

5(2s, w):= ^ exp[u;C(n, v)] = C{2s) F,,^ o (/ - H,,^)"! [1] (0) . 

(M,D)en 

Using well-known properties of the Riemann zeta function Ci^): all our re- 
sults for Oat will follow from those on f^Tv- 

In view of (EH ElSl EHEl) , the relations (|TT7j) . and (ITTHl) connecting the 
Dirichlet moment generating function with the transfer operator are the 
analogues for rational trajectories of the relation (|2.1()|) for the truncated 
real trajectories. In the case of rational trajectories, we have to work with 
the quasi-inverse and extract the coefficients of Dirichlet series: This is why 
the discrete problem is more difficult to solve than the continuous problem. 

2.5. Perron's formula and Dolgopyat's estimates. We wish to evaluate 
the sum ^^iN) of the first coefficients of the Dirichlet series S{2s,w). 
Our first main tool towards this goal is the Perron formula. The Perron 
formula of order two (see e.g. jl?!) applied term by term to a Dirichlet 
series F{s) = J2n>i ^nn~^ and a vertical line = D > inside the domain 
of convergence of F says that 

(2.19) ^(T) := ^ a„(r - n) = — F{s)——- ds . 
Applying Perron's formula to the Dirichlet series S{2s,w), we find 

1 i-D+ioo rp2s+l 

(2.20) m^{T):=Y,Cn{w){T-n) = — S{2s,w) ds . 

Thus, Perron's formula gives us information on ^'^^.(A^), which is just a 
Cesaro sum of the ^wiQ)- 

*^(iv) =Y.Y1 = J2 ^-(^) ■ 

Q<Nn<Q Q<N 

Proposition and Relation 1)2. 17() show that s i-^ S{2s, w) has a possible pole 
at s = cr(w), where a{w) is the unique complex number near 1 for which 
X{a{w),w) = 1. In particular, the integral in (|2.2U|) is well-defined only 
if D > sup^(5?(j('t/;)) (note that supy^{^a{w)) > 1 can be made arbitrarily 
close to 1 by taking w close enough to zero). To combine the Perron formula 
with Cauchy's residue theorem, we wish to modify the integration contour 
= D into a contour containing a{w) as a unique pole of S{2s,w). This 
is possible and leads to a quasi-power expansion if there is a > 0, such that 
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(i) S{2s,w) admits s = a{'w) as a unique pole in the strip \?R.s — l\ <a. 
(a) In a "perforated" half-plane {?R.s > 1 — oq , |s — 1| > ao/2}, for small 
a > ao > 0, the estimates \S{2s,w)\ < max(l,(9s)^) for < ^ < 1, hold, 
uniformly in w close to 0. 

Note that (i) cannot be satisfied if the map T is conjugated with a 
piecewise affine map, since, in this case the Dirichlet series S{s, 0) has an 
infinite number of poles in the vertical 'Sis = 1. Note also that bounds of 
the type (m) are extremely difficult to obtain for general Dirichlet series. By 
the above discussion, such bounds would follow from similar estimates on 
the quasi-inverse of the operator, which are closely related to those obtained 
by Dolgopyat ^Hl- In the spirit of Chernov [2,, Dolgopyat 16 introduced 
several "uniform nonintegrability" ( UNT) conditions. They allowed him to 
control oscillatory integrals associated to iterates of transfer operators for 
s = a+it, with fixed a close to 1, and prove exponential decay of correlations 
for some flows. We shall give a new formulation of Dolgopyat 's strongest 
such ?7A/^/ condition, which appeared implicitly in Section 5 of his paper |16j 
and turns out to be satisfied by our three algorithms. This condition is 
stated as an assumption on the derivatives of the inverse branches of the 
dynamical system in Section 3.2, and expresses that, in a sense, the map T 
is quite different from a piecewise affine map. In Section 3, we shall prove 
the following theorem, which is the central functional analytic result of the 
paper. In the statement, we use the following family of equivalent norms on 
CHl): 

(2.21) ||/||i,t:=sup|/| + ^^^, t^O. 

Theorem 2 [Dolgopyat-type estimates]. Let {I,T,c) he a triple ofQMQ- 
type, with contraction ratio p < 1, and such that Condition UNI from §3.2 
holds. Let iis,w be its weighted transfer operator i|2.6'|) acting on C^(T). 
For any ^, with < ^ < 1/5, there is a (real) neighborhood Si =]1 — a, H-a[ 
of 1 (which depends only on (T, T) and not on c), a (real) neighborhood Wi 
of 0, and there is M > such that, for all s = a + it, w = v + it with 
(a,iy) G Si X VFi and \t\ > l/p^ 

(2.22) ||(/-H,,^)-ii,t < M- 

Note that we shall have to modify Dolgopyat's arguments since we must 
consider dynamical systems which possess an infinite number of branches 
(see in particular Lemma 1), and we work with bivariate weighted transfer 
operators H^^^ involving a cost function. 

2.6. Statement of the Central and Local Limit Theorems. We shall 
see in §4.1 that Perron's Formula 1)2. 2U() combined with the fundamental 
relation 1)2. 17() . together with the bounds a la Dolgopyat (Theorem 2) pro- 
vide a quasi-powers estimate for the Cesaro sum ^'^„(A^). It does not seem 
easy to transfer this information on ^yi;{N) to estimates on ^yi]{N), because 
the coefficients are complex. The way we overcome this is by first proving 
(Lemma 11 in Section 4.2) quasi-power estimates for the moment generating 
function of some "smoothed" version of the cost C, for which the transfer 
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is possible by standard methods. We are then able to apply Theorem 
to the smoothed model, and show that the two models are close enough in 
distribution so that the following holds: 

Theorem 3. [Central Limit Theorem for rational trajectories.] For a Eu- 
clidean algorithm amongst Q, JC, O, there is 7 > 0, so that, for any cost c 
of moderate growth, letting A(s) he the function from Proposition 0: 
(a) The distribution of the total cost C on Vlj^ is asymptotically Gauss- 
ian, with speed of convergence 0(l/\/Tog7V), i.e., there exist two constants 
fj,{c) > and 5{c) > such that, for any N, and any y G M 



N 



I X C{u,v) - /i(c) log iV 
o(c)vlogiV 



1 o , /I 



(6) The mean and the variance satisfy E]\[[C] = fi{c) log N + ry(c) + 0{N~'^'), 
and Yn[C] = 6^ (c) log N + 6i{c) + 0{N-^). 

Generally, for each k > 1, there is a polynomial Pk of degree k so that 

E^[C']=P,{logN) + o(^^^^' 
with a O-term uniform in k. 

(c) In the special case c= I, denoting fi := /i(l), 6^ := 5^(1), we have 

2 .2 2|A"(1)| 

In the general case, 

fi{c) = fi ■ Jl{c) , 5\c) = /22(c) -d^ + fi- ?(c) + fi^Jl{c) ■ x{c) > , 

where /1(c) > and S^{c) > are given in Theorem 1, and x(c) = A"^(l, 0). 
Claims (a), (b), and (c) also hold for Pjv on ri^r. 

Note that 2/fi is the Kolmogorov entropy for (T,fidx). Also, 7 does not 
depend on the cost. The constant x(c) can be viewed as a covariance coef- 
ficient between the number of steps P and the cost c. Since there exists a 
closed form for fi in the three cases of interest (cf. Figure 1), the constants 
jj,, and thus /i(c) can be easily computed (see remark after Theorem 1 and 
j49j ) . The constants 5, 5{c) are proven to be non zero in Proposition 1, 
Section 3.5. They do not seem to admit a closed form. However, Lhote has 
proved that they can be computed in polynomial time [32j . 

In Section 5, restricting to lattice costs, we obtain bounds for K]\f[exp{iTC)] 
with T G [— vr, vr], and combine them with saddle-point estimates to get: 

Theorem 4. [Local Limit Theorem for lattice costs.] For any algorithm 
among Q, K., O, and any lattice cost c of span L and of moderate growth, let- 
ting 11(c) > and 5'^{c) > be the constants from Theorem 3, the following 
holds, 



N 



(u, v); < C{u, v) - /i(c) log N - (5(c)x0ogiV < ^ 
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(5(c) V27r log iV \logNj' 
with a O uniform for x E M. The same holds for Pjv in ^n- 




3. Property UNI and Dolgopyat estimates. 



In this section we shall first prove Theorem 2 and check (§3.5) that its ad- 
ditional assumption (UNI) holds for our three Euclidean algorithms. Then, 
we show that Theorem 2 together with aperiodicity results imply useful 
estimates on S{s,w). 

To prove Theorem 2, we use ideas due to Dolgopyat ^Hl; however, we have to 
adapt them to our context, i.e., bivariate weighted transfer operators H^ -^; 
associated to triples (I,T,c) of ^A^^-type with infinitely many branches, 
as explained after the statement of Theorem 2 in Section 2. One of the main 
ideas of Dolgopyat was to deal first with the L^-norm of some iterate H"^ [/] , 
with an index n that depends on t = 3=5. Then, he made two transfers of 
estimates: first from this L^-bound into a bound for the sup-norm, next from 
this sup-norm-bound into the desired bound for the (l,t)-norm. Following 
his strategy, we establish preliminary results in Lemmata 1, 2 and 3. (Note 
that Lemma 1 is new.) §3.2 is central: our version of the UNI Condition 
is stated and shown to entail the desired exponential estimate for the L^- 
norm of the operator (Lemmata 4 and 5). The two transfers are operated 
in §3.3: they lead to Theorem 2. We then check in Lemma 6 (§3.4) that 
its additional assumption (UNI) holds for our three Euclidean algorithms. 
In §3.5, we show both an aperiodicity result on vertical s-lines and ri;-lines 
together with a convexity result (Lemma 7 and Proposition 1). Finally, in 
§3.6, we return to the Dirichlet series S{s, w) and obtain further estimates in 
the compact neighborhood of the real axis which is not covered by Theorem 
2 (Lemmata 8 and 9). 

3.1. Preparatory material. Transfer operators iis,w and their nor- 
malizations iis^w Triples (T, T, c) of QMG type and their associated 
transfer operators ils,w, acting on C^(X), for {^s,^w) G Sq x Wq were 
introduced in § 2.2. We summarized in Proposition well-known spectral 
properties that we shall need, in particular the existence of a dominant eigen- 
value A((T, z^), [with positive eigenfunction fa-^u] or X{s,w) for suitable real 
{(7,1') G So X Wq [from 1)2. 5|) ] and complex (sjw) close to (1,0). Recall that 
is the set of inverse branches of T"', TC* the set of all inverse branches of 
any depth, p < 1 the contraction constant and K the distorsion constant. 

Defining K = K/{1 — p), it is not difficult to check that. 



The above bounded distortion property will play an important role. Note for 
further use that, since the ratio h"{x)/h'{x) is the derivative of log|/i'(x)|, 
setting L := e^, we get 



h" {x)\ < K\h' {x)\ , yx £l,h£n* . 



(3.1) 



1 \h'{x)\ 

L-WM 



< L 



for all x,y I, for all h £ TC* . 
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It will be convenient to work with the normalized operators iis^w defined by 

(3.2) U.s ^[f] = — — 3— — Ils^w[fa,u ■ f] , s = a + it, w = u + ir . 

By construction, for {(T,^) G Sq x Wq, the operator ila,u acting on C^{T) 
has spectral radius equal to 1, and fixes the constant function = 1. Also, 
Sp'H.a+it,u = ^{o',i')Sp'H.a+it,u- Remark next the inequality ||Hs,^[/]||o < 
|oHo-,;/[l] = II/IIoj which implies the useful bound 



(3.3) ||Hs^tt,||o < 1 . 

It is easy to check that H*_^ fixes the probability measure fia,u = fa,u • V-cr,u- 
Remark about notations. In the sequel, the notation A{x) « B{x) 
means: A is less than B up to absolute multiplicative constants. This means 
that there exists some absolute constant k such that for every x of interest, 
A(x) < kB{x). It is synonymous with A(x) = 0{B[x)) with an absolute O- 
term. The symbol W denotes a complex neighborhood of for the variable 
w. If J Cl is a union of intervals, we denote by \J\ its Lebesgue measure. 

Relating Ho-,^ and Hi q. In order to exploit properties of Lebesgue mea- 
sure, which is fixed (only) by the dual of ilifl, Dolgopyat uses the following 
property (see e.g. last lines of p. 367 in ^^). When (I, T) has finitely many 
branches, there is A^ — > 1 as o" — > 1 so that for positive / E C^(X) 

(3.4) H,,o[/](^) <^aHi,o[/](x). 

The above inequality is not true in general when there are infinitely many 
branches (it fails for the Gauss map). The purpose of the following lemma 
is to extend (jSHJ to the case of infinitely many branches and bivariate 
operators, comparing H"j^, and to their analogues for {(t,v) = (1,0): 

Lemma 1. For {a,v) € So x Wq, denote 

__ A(2a- 1,2^)1/2 

Let C be a compact subset of Sq x Wq. For J C , denote by J = 
Uhejh{I). Then, for (cr, z^) G C, 

furthermore, for any f £ C^{Z), for any integer k > 1, 

(3.5) ||H^,,[/]||2«a2^||hJ^o[|/|2]||o. 
The absolute constants involved only depend on L. 

The function A^j^y depends continuously on {a,v) and A\i^ = 1. 

Proof. The equality Ha^vlf] = /^cr,!/[H^^[/]], when applied to the character- 
istic function of some fundamental interval h{Z) of depth k proves that 

ii^Am\ « ^^^T^ / \h'{x)rdf,^A^). 
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Moreover, by the bounded distortion property p.lj) . the ratios (two by two) 
of the three quantities a(/i), c(/i), 

a{h):= j |/i'(x)r(i^^,^(x); := |/i(J)r ; c(/i) := j \h' {x)\'' dx , 

admit upper and lower bounds that do not depend on /i, and are uniform 
for {<J,v) € C. Then, summing the inequahties 

exp[i/c(/i)] 
\{a,v) 

over J C Ti^ ^ and applying the Cauchy-Schwarz inequality, one gets 

^i.AJ\ < E exp[2z.c(M] • \h{i)\''^-'^ \h{i)\^ . 

Then, dominant spectral properties, together with bounded distortion, entail 
the inequality 

E exp[2uc{h)] ■ \h{I)\'^''-^ « X{2a - \,2vf , 

and, finally, the relation | J| = J2heJ provides the first claim. 

Consider now / G C^iT). The relation 

l(H.,.)^[/](x)|«— ^ E expK(/i)]-|/i'(x)r-|/oMx)|, 
is valid if (fi, v) belongs to £, and, by the Cauchy-Schwarz inequality. 



2 



E expK(/i)]-|/i'(x)r-|/o/i(x)| 



The second factor is exactly HJq[|/P](x), which is less than H^' Q[|/p](a:;) 
(up to absolute multiplicative constants). Thanks to dominant spectral 
properties, the first factor is easily related to \{2a — 1, 2v)^ . □ 

Lasota— Yorke bounds. The following lemma describes how Ylg^w acts 
with respect to the quasi- norm ||.||i when s varies over a vertical line: 

Lemma 2. For every compact subset £ of Sq x Wq, there is C > 0, so that 
for all {s,w) with {?fis,?fiw) £ C, and all f £ C^{I), 

lfc/||i<C(|.|||/||o + p"||/||i) , Vn>l. 



Proof. The quantity H" ^[/] can be written as a sum over h £ TC^ of terms 

exp[ifc(/i)l , , . , ,, 1 , ^ r 

w J rh{x) with rh := • ^ • iUuf) o h . 
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The Leibniz sum for the derivative of contains three terms. We can bound 
the first for all s using the distortion assumption since 

< \s\K\\h'\'\ = |s|i^|/iT- 

Compactness of C and continuity of (a, i^) i— > fa-^u, and (cr, i/) ^ f'^y imply 
that the second term may be controlled by 

If I 1 



fa,w fo;w 

for some Cc > 0. Finally the last term can be estimated using 

KU^ ■ /)' o h\\h'\ < ■f\oh + iu,^ ■ i/'i) o h] . 

We can ensure 

K\s\ + Cc + p''Cc<C{C)\s\, 
so that the derivative (H"^[/])' satisfies 

||(H^,^[/])'||o <c{\s\ ||H^,J/]||o + pn\U:A\f'\]\\o 
The final result follows from □ 

First use of the (l,t)-norm. In the bound from Lemma 2 for the 
derivative of H"^/(x), there appear two terms, one which contains a factor 
|s|, the other a decreasing exponential in n. In order to suppress the effect 
of the factor \s\, Dolgopyat uses the family of equivalent norms 

II/IIm ■= ll/llo + ^ll/lli = sup|/| + lsup|/'|,t / 0, 

which appear in the statement of Theorem 2. With this norm and Lemma 
2, together with (|3.3|) . we obtain the first (easy) result: 

Lemma 3. For any ti > 0, for every compact subset C of Sq x Wq, there is 
Mq > so that for aU n> 1, all {s, w) for which "Siw) E C and \ Qs\ > ti 
we have ||H^^^||i,q^ < Mq. 

3.2. UNI Condition and L^-estimates. Assuming UNI, Dolgopyat first 
proves that there is 7 < 1 so that 

|Hr[/](^)l'^x<7"1l/llM, 

for all large t and hq = 0(log \t\). In this subsection, we extend this result to 
the bivariate operator H^^^, when the number of branches of T is possibly 
infinite. 

Writing s = a + it, w = u + it, |H"„,[/](x)p can be expressed as: 

(3-6) ^ -.orj exp[wc{h) +wc{k)] ■ exp[it^'/, ^(x)] • Rh,k{x) , 

^ ' (/i,fc)6'H"x-H" 

with 



(3.7) ^ft,fc(x) := log 



\k'{x)\ ' 
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(3.8) RhA^) = \h'{xr\k'{xr-^—{f ■ /,,,) o h{x) ■ if ■ /,,,) o k{x) . 

Using / = 3f?/ + iQf, the term Rh^k decomposes into four terms, each of 
which has the form 

(3.9) rh,k{x)=e'^\h'ixr\k'ix)\''-^ig-U^,)ohix)ii-U,,)ok{x), 
for two real functions g,i £ 9/} and expiw G 

The functions ^h,k play an important role here: the sum 1)3. 6() will be split 
into two parts, according to their properties. The first sum will gather 
the pairs for which the derivative |^/ifc(a;)| has a "small" lower bound, and 
condition UNI will precisely require that there are not "too many" such 
pairs {h, k), providing a convenient bound for the corresponding integral /~ 
(Lemma 4). The second sum will gather the other pairs {h, k), for which the 
derivative |^'^^(x)| has a "large" lower bound. In this case, the Van Der 
Corput Lemma on oscillatory integrals will be applicable (Lemma 5), giving 
a bound for the corresponding integral I A 

Let us introduce some notations needed for our formulation of UNI. For two 
inverse branches /i et /c of same depth, we introduce a "distance:" 



(3.10) A(/^,fc) = inf |^^(x)| = inf 



For h in TL"', and t] > 0, we denote 

(3.11) Jih,v)-= U ^(^)- 

keH" ,A{h,k)<ri 

Property UNI(a) expresses that the Lebesgue measure of J{h,S) is « S 
when 5 is scaled similarly to the maximal length of fundamental intervals 
of depth n. For any p > p, this length is 0{p^) (up to absolute constants) 
and plays the role of a reference scale. This is a reformulation of the UNI 
condition implicit in Dolgopyat's Section 5 [16,, which we (finally) state: 

UNI Condition. A dynamical system of the good class, with contraction 
ratio p, fulfills the UNI Condition if each inverse branch of T extends to a 
function and 

(a) For any a (0 < a <1) we have | J(/i, /9'='")| « p""" .Mn^Mh £ m . 

(b) Q :=sup{|^'';^ ;-(x)|;n > 1 , /i, € x e J} < oo. 

Remarks. Note first that UNI does not involve the cost: this is because 
c is constant on the monotonicity intervals of T. 

For dynamical systems with affine branches, all the A(/i, k) are zero, and, 
for any r/ > and any h G 7i*, the interval J{h, r]) equals 2. Then dynamical 
systems with affine branches cannot satisfy UNI. We will see in Proposition 
1, Section 3.5, that this is the same when the map T is conjugated with a 
piecewise affine map. 

Condition (b) follows from the existence of Q < cxd so that 

(3.12) \h"'{x)\ < Q\h'{x)\ , Vn > 1 , V/i e . 
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It suffices to check ()3.1'2() for n = 1 (similarly as for the distortion condition) . 
(Note that this condition is always satisfied if there are finitely many 
inverse branches). 

Study of the L^-norm: the close pairs. 

Lemma 4. Recall A^^y from Lemma 1. Suppose that Condition UNI(a) 
holds. For any compact subset £ of Eq x Wq, for all {a = 5Rs, u = Rw) £ C, 
for all n, for all a, with < a < 1, the integral I~ of the sum p.6'|) restricted 
to pairs {h, k) G Ti^ x Ti^ for which A(/i, k) < satisfies 

\I-\ = \I-{s,w,f,a)\ « 



Proof. Up to a positive constant that only depends on (a, v) (and is uniform 
on the compact subset £), is less than 



211 



exp[u{c{h)+c{k))]- f \h' {x)\'' \k' {x)\'' dx . 



First, using the bounded distortion property, for all pairs {h, k) G Ti.* x Ti.* 
up to multiplicative absolute constants, one has 



h'{x)\''\k'{x)\''dx « i^j \h'{x)\''dx 



\k'{x)\''dx 



Then, as in the beginning of the proof of Lemma 1, with the same bounded 
distortion property and using that is an invariant probability for the 
normalized operator, the ratios (two by two) of the four quantities a(/i), 
b{h), c{h), d{h), 

'^{(^t'^) Ji '^[c,'^) Ji 

c{h) := f^^Am] ; d{h) := ^^EMM|/,(T)r ; 

A((j, v)^ 

admit upper and lower bounds which do not depend on h and are uniform 
when (cJ, i^) varies in a compact set. Up to a multiplicative constant, it is 
then sufficient to study the sum 



heH" 



\ 



ken" 



heH" 



\ A{h,k)<p"-"- 

Now, the first relation of Lemma 1, fia,u[J] ^ C^^u^^u l-^l^^^' which holds for 
any subset J that is a union of fundamental intervals of depth n, is applied 
to J{h,p"''^). UNI (a) provides an evaluation of its Lebesgue measure, and, 
finally, |/-| « (p^/^A,,,,)" 1 1/| Ig. □ 



Study of the L^-norm: Application of the Van der Corput Lemma. 

Consider now the integral of the sum relative to pairs (/i. A:) which were 
not treated by Lemma 4: 
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Lemma 5. Suppose that Condition \JNl(b) holds. Letting [x] denote the 
smallest integer > x, set 

1 



(3.13) 



no = no{t) 



log p\ 



log |t| 



Then, for any compact subset £ of Sq x Wq, for any {a = "Sis, v = "Siw) G £, 
and \t\ = \Qs\ > 1/p^, for any < a < 1/2, the integral of the sum l^3.6\) 
for n = no, restricted to {h, k) £ H"" x with A{h, k) > /^''"o, satisfies 



2 

l.t ■ 



Proof. We start with a general integer n and bound 



J] exp[zy(c(/i) + c(A:))] • I{h, k) 



A{h,k)>p"-"' 

where the integral I{h,k) := jjex.p[it'^h,k{x)] Rh,k{x)dx involves '^h,k, Rh,k 
defined in (|3.7|) . (|3.8j) . and decomposes into four integrals of the form 

I{h,k) := j exp[it'I^h,k{^)]^h,kix)dx , 

with rh^k defined in (|3.9)1 . We shall apply the following lemma to each 
oscillatory integral I(h,k): 

Van der Corput Lemma (See e.g. |lSj). For each interval 2 and every 
Q >0, there is C{Q), so that for alJ t e M, ^' E C'^{1) with {x)\ < Q , 
|^''(x)| > A with < A < 1, and r £ C^{T) with ||r||o < -R, ||r||i,i < 
RD, the integral I{t) = exp[it^'(j;)] r{x) dx satisfies 



\Iit)\<RC{Q) 



D + l 1 



Consider (t, n) with l/\t\ < p''". Setting 



1 



(3.15) M(/i,A:) :=sup|/i'(x)nA:'(x)r ^„ , , 
xex J a Ax) 



fa,u O h{x)f^^y O k{x) , 



the norm ||r/i^fc||o satisfies 

||r/,,fc||o <M(/i,A:)||g||o m\o < M{h,k)\\9\\i,t IKHm- 

The arguments used in the proof of Lemma 2 for the function apply to 
the function r/j^^, and 

\\rhA\i,i « Mih, k) [Ibllo (ll^llo + IKIli) + ll^llo (llffllo + P" WoWi)] 

« M(/i,A;)||<7||M||^lli,t[l+P"|t|] • 

Then, by Property UNI(h), the Van Der Corput Lemma can be applied to 
each integral I{h,k), which thus satisfies 



\Iih,k)\ «M{h,k) |bl|i,t||^||i,t 



2+t/3'" 1 

, , + 



,2an 
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Now, we choose n = no as in H3.13() . Since a < 1/2, and > 1/p we have 
no > 2 and < p~("o-i) < |^|^ gQ that the Van der Corput Lemma may 

be apphed. The previous inequahty becomes 

(3.16) \I{h,k)\ «M{h,k) ||<7||m||^||m/9(i-2'^)"«. 

Returning to the integral I{h, k), 

\T{h,k)\ « M{h,k) p{l"2a)no 11^112 



l,t ■ 

Now, take xq in I. Then, from the bounded distortion property 1)3. Ij) . and 
the definition of M{h, k) in (jSHSJ), we get 

M(/i, k) « \h'{xor\k'{xo)r ^ ] U,u o h{xo)U, o k{xo) , 

and therefore 
(3.17) 

— ^ ^ eMHc{h) + c{k))]M{h,k) « (h^,,[1](xo))' = 1. 



From (ITTfill3J71l3T^ . we finally obtain |/+ | « p(i-2")'^o \\f\\l^. □ 

Study of the L^-norm: the final result. Consider the integer uq from 
(|TT^ of Lemma 5 (for \t\ > and some a with (2/5) <a< (1/2). Then, 

since a/2 > (1 — 2a) > 0, there exists a (real) neighborhood of {a, v) = (1, 0) 
on which 

(3.18) A^^u ■ /o"/^ < P^"^" for any (a, v)eW x^. 

Then, from Lemmata 4 and 5, 



(3.19) j |H^;^[/](x)pdx «p(i-2«)"o 



2 



3.3. End of proof of Theorem 2. We operate now the transfers between 



various norms. 

2 



From the L -norm to the sup-norm. Since the normalized density 
transformer Hi is quasi-compact with respect to the (1, l)-norm, and fixes 
the constant function 1, it satisfies 

(3.20) \\&i[\g\^]\\^=(^j^\g\\x)dx^ + OirDWg'^W^,^ , 

where ri is the subdominant spectral radius of Hi. 
Consider an iterate H" ^ with n > uq. Then 

muml « iiH^;""b]iig with g = lu^ium . 

Now, using (|3.5j) from Lemma 1 and (|3.20j) with k := n — no, together with 
the bound (|3.19j) for the L^-norm and finally Lemma 2 to evaluate Hff^Hi^i, 
one obtains 



p(l-2a)no _|_ ^n-no |^| 



2 

l,t 
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We now choose n = ni as a function of t so that the two terms p(i-2a)no 
and r"""""!*! are almost equal (with no(t) defined in H3.13|) l: 

(3.21) ni = (l + 7])no with r? := 2(1 - a)-^^^ > . 

logri 

Choose now d such that < 7y(5a — 2) <d<l — 2a<l/5 (which is possible 
if a is of the form a = 2/5 + e, with a small e > 0). Recalling 1)3.18(1 where a 
first neighborhood was defined, and considering a (real) neighborhood T,xW 
of (1,0) for which 

(3.22) sup [A(cT, z.)i+'', Al^] < p-'?(W2-i) < p-d/2 ^ 



we finally obtain, for ni(t) and rj defined in ()3.2H) 

(3.23) ||H?,1.[/]||o«p"^'||/||m, with 6:=i^^^ 



From the sup- norm to the ||.||i j-norm. Applying Lemma 2 twice and 
using H3.23() yields the inequality 

r2n 



Hf^[/]||i « |.|||H"^^[/]||o + />"M|H"i„[/]||i 



« |.|p"i^||/||i,, + p«i|t|(^^||/||o + p"i 

(3.24) « Itlp^^^ll/lli.t, 

which finally entails for n2 = 2ni (and ni(t) as above) 

(3.25) ||H^,L||i,t«p"^^/2. 

The last step in Theorem 2. For fixed t with |t| > 1/p^, any integer n 
can be written n = /cn2 + i with £ < n2{t). Then (|3.25|) and Lemma 3 entail 

\\^lJ\i,t < Mo llH^f^lltt < Mop^'''''/^ < Mop^"/2^-^n2/2 _ 

Since bn2/2 = bni = (1 — 2a — d)nQ, with no defined in ()3.13|1 . we finally 
obtain 

\\nlJ\^,t<M\t\^^\ 

with ^ := 1 -2a-d, 6:=^—, 7 := p''/^ < 1, Mi = 

1 + 77 

Then ^ is any value between and 1/5. Therefore, returning to the operator 
^, we have shown 

(3.26) ||H" ||i,t<Mi-7"-|t|«-A(cT,i/)", Vn , V|t| > 1//?^ . 



Finally, for any (a, i^) G S x 1^ as in (|3.22() . one has 

i i i 

7A(CJ, I^) < + • p 4(l + r,) = p1(l + v) = 7 < 1 . 

This proves Theorem 2 with M := Mi/(1 — 7). □ 
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3.4. UNI Condition and Euclidean dynamical systems. To apply 
Theorem 2 to our three algorithms, we prove that they satisfy UNI. 

For two LFT's hi and /i2, with /ij(x) = (a^x + bi)/{ciX + di), we have 

\c1d2 - C2di\ 



h" 
h' 



x) 



h" 



|(cix + di){c2X + d2) 



so that the distortion property gives 
(3.27) A(/ii,/i2) 



di 



C2_ 
d2 



inf 



h[{x)h'^{x) 



1/2 



> 



I 



di 



C2 
(^2 



Hence A only depends on the difference of the quotients Ci/di of the denom- 
inators of the LFT's. We shall next show that this difference is the "honest" 
(ordinary) distance between the rationals /i^(0) and /i2(0), where h* is the 
mirror IFT of h, defined by: 



h*{x) 



ax + c 



if 



h{x) 



ax + b 



bx + d ^ ' cx + d 

This mirror operation appears in |42j where Schweiger relates it to the nat- 
ural extension, and in [7j, where the authors use the geometric notion of 
"folded" and "unfolded." 

Clearly, the mirror map is an involution satisfying the morphism prop- 
erty (hok)* = k* o h*. It is not difficult to see that ii h £ TiP 
is a LFT from a Euclidean dynamical system, associated to the se- 
quence eo, (mi, ei), (m2, €2) ■ ■ ■ (m-p, e^), the mirror LFT h* corresponds to 



[m 



p— 1 ) ep_2) . . . (mi, eo), i.e., the decomposition involves the 



same "digits" as h, but in the inverse order. 

By ()3.27|) . the distance A(/ii, /12) between hi and /12 indeed satisfies 



(3.28) 



A(/ii,^2) >Y\hm-him- 



It is not difficult to check for each of our three Euclidean algorithms that 
the set {h*,h £ 7i} is the set of inverse branches of a dynamical system 
(T*,r*), the dual dynamical system, which belongs to the good class, with 
the same contraction ratio p* = p and a distortion constant I*, as we now 
explain: For the Classical Euclidean algorithm Q, since all the e are equal 
to 1 note that (I*,T*) = (2,T). For the two others, the equality p = p* 
follows from the definition (|2.4|1 of p, the distortion property for the dual 
system, and the fact that, in both systems the worse branch, for which the 
contraction ratio is attained on a fixed point, is a LFT h with h = h* . The 
three dual dynamical systems are described in Figure 2. 

Our next goal is to check UNI for our three algorithms. Note first that 
checking UNI (b) amounts to verifying p.l2p . which is easy for the three 
algorithms. For UNI (a): 

Lemma 6. A Euclidean dynamical system of the good class which admits 
a dual system that belongs to the good class satisfies UNI (a). In particular, 
the dynamical systems associated to the Euclidean algorithms Q,JC, and O 
satisfy UNI. 
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Algorithm 


g* 


fC* 


o* 


Interval 


[0,1] 


[cb- 


2,0-1] 






-2,0 




Generic 


TO > 1 ,e = +1 




m > 2,e = ±1 


TO > 1 odd, e = ±1 


conditions 








if TO 


= 2 then e = +1 


if m = 1 then e = +1 


































1 !1 \ 
)\\\ \ 








i 

\ 




Graphs 





















Figure 2. The three dual Euchdean dynamical systems: 
Standard, Centered, Odd. (0 = 



Proof. Fix < a < 1 and h E H'^. Denote by J*{h,r]) the union of the 
intervals k*(I*) for LFT's k G Ti.'^ satisfying A(/i, k) < r]. First, we estimate 
\J*{h,r])\, second, we relate \J*{h,r])\ and \J{h,r])\. 

By (HT^l) if A{h,k) < r] then \h*{0) - k*{0)\ < Lrj. Since the length of a 
fundamental interval of depth n in the dual system is at most (up to an 
absolute constant) we get 

(3.29) |J*(/i,r/)l « 2L?7 + 2p''". 

The fundamental intervals k{I) and k*{I*) (of depth n) have almost 
the same length: Indeed, by the bounded distortion property, and since 
(r)'(O) = A:'(0) for all k e n\ 

^ ' (LL*)i/2 - ]fc*(T*)l lfc'(0)| lfc*(T*)l - ^ ^ • 

Thus, since the intervals k*{Z*) in J*{h,rj) are disjoint, if r/ < p"", using 

(3.31) |J(/i,?7)| < (LL*)i/V*(^^)l « (LL*)^/^(2L?7 + p«") « p-^" . □ 

3.5. Condition UNI, aperiodicity and absolute convexity. Condition 
t/7V7 provides estimates for the (l,t)-norm of (/ — H^^io)^^ when \t\ = \Qs\ 
is sufficiently large. We have now to consider the case when jtj is not large, 
and we explain how the f/iV/ Condition intervenes in this context, via aperi- 
odicity results. We must also check that variance constants 6{c) which will 
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appear in Theorems 3 and 4 are not zero, and the UNI Condition intervenes 
in this context too, via absolute convexity results. 
We first recall a classical result. 

Lemma 7. [Aperiodicity and absolute convexity.] Let Ylg^w be the transfer 
operator associated to {I, T, c) of QM.Q type. 

(i) [Aperiodicity.] Denote by R{s,w) its spectral radius, and consider a 
point {to, To) 7^ (0, 0). The following are equivalent: 
(a) There exists (cr, i/) G Sq x Wq for which 

R{<T,u) e SpHo-+it(,,,y+iro • 

(6) There exists f G C^{I), \ f\ = 1, so that, for aU n, for all h^W 

(3.32) |/i'(x)P*o •exp[iroc(/i)] •/o/i(a;) =/(x), VxGT. 

{ii) [Absolute convexity at (1,0)]. Consider a point {q,r) ^ (0,0), the op- 
erator Hi+g^^rw and its pressure function w i-^ A(l + qw, rw) . The following 
are equivalent 

(a) The second derivative of w ^ A(l + qw, rw) is zero at w = 
(6) There exists a function f E C^i^) strictly positive and a constant 
a > , so that, for all n, for all h G TC" 

(3.33) |/i'(a;)|9-exp[rc(/i)] •/o/i(x) =a"-/(x), VxGT. 

Proof, (i) See for instance Prop. 6.2 in or Prop. 9 in 45 . Since 
Sp Ho-_|_jtQ^jy+jT- = X{a,iy)Sp'H.a+it,u we may replace R{a + it,!/ + it) and 
i?((T, v) by the spectral radii of the corresponding normalized operators. By 
Proposition 0, the spectral radius of 'H.a+it,u+iT is at most 1 = R{H.a,u)-, 
while its essential spectral radius is at most p < 1. If (a) holds, then 
'^(T+ito,u+iTo has an eigenvalue A = 1 with an eigenfunction / G ^^{1) with 
max/ I /I = 1. Suppose that this maximum is attained at xq G X. Then the 
equality W^j^^f.^^^j^^^^^[f]{xQ) = /(xq) can be written as Y^haH^ ^hhh = 1 with 

ah ■= -r—? — T exp[i/c(/i)] • f^^^ o h{xo) ■ |/i'(xo)r 
^ -/o/i(xo)exp[iroc(/i)]|/i'(xo)r*° 



Normalization implies that Y^heH" = ^ while each factor has modulus 
at most 1. Therefore, each equals 1, and, for every n > 1, 

f o h{xo) = /(xo) exp[-iroc(/i)]|/i'(xo)r'*° , V/i G . 

In particular, by density, |/| is the constant function 1 and the above remarks 
also hold for any x G X. Then, for all x G /, all n > 1, and all h ^ TC, 

(3.34) jM_ = exp[.roc(/i)]|/i'(x)r*o . 

Following the arguments backwards, we get the other implication. 
(a) We apply to the operator Hi+g^^^to the results of [HI. In Proposition 6.1, 
Broise states her results in the context of functions with bounded variation. 
Due to our strong Markov assumption, we may work in C^(X), and the 
cost of interest is D := — glog|r'| + rc. Since fi is a strictly positive 
function and Hi[Z)] belongs to C^, we may transfer Broise's proof to our 
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context : it shows that Condition {ii){a) is equivalent to the following: 
there exist u G and a constant K for which, for any h G 7i, one has 

glog \h'\ + rc{h) = u — uoh + K. This last condition is clearly equivalent 
to {ii)(b) with / := exp[n]. Note that the fonction u introduced by Broise 
which involves the centered version D of D, i.e., D := D — ijl{D) 

u:=-1(7-Hi)-1oHi[;D-/i] 
/i 

actually belongs to C^(X). 

□ 

Conditions {i){h) and {ii){h) are of the same form and closely related to con- 
jugaison with piecewise affine maps. As we next see it, the UNI condition 
entails that this conjugaison cannot occur. Then, under the f/A'^/ condition, 
the variance constants will be always strictly positive. On the other hand, 
the aperiodicity result, together with the LW/ condition, provides some use- 
ful informations about the spectrum of Hs^„,., notably for lattice costs. We 
recall the definition of lattice costs: a cost c is said to be lattice if it is not 
zero, and there exists L > for which c/L is integer. The largest such L is 
the span of the cost. 

Proposition 1. Consider a QM.Q system satisfying UNI(a). The following 
holds: 

(i) The mapping T is not conjugated with a piecewise affine map. 

(ii) The pressure function of the operator Jli^qu,,rw is absolutely convex at 
w = for all fixed r G M and q ^ 0, (i.e., its second derivative is strictly 
positive). In particular, for r = 0, one has A"(l) > 0. 

{in) For any t ^ 0, 1 does not belong to SpHi+j^^o- Furthermore, if c is 
lattice with span L, for all t, and all r not multiple of 21: /L, 1 does not 
belong to SpHi+u^r- 

Proof, (i). Suppose that T is conjugated with a piecewise afHne map. 
Then, there exists / > in C^iT) such that, for any n, for each h G H", 
there is a constant d{h) for which \h' {x)\f o h{x) = d{h)f[x) for any x G T. 
Then taking the logarithm, differentiating and putting / := log / G 
we get 

If ;^^(,t) = [h'{x)f' o h{x) - k'{x)f' o k{x)] , Vn , V/i , yfc G . 

Then, for p < < 1, the distance A(/i, k) satisfies A(/i, k) « p" for any h, 
k eW^, which contradicts UNI (a). 

(ii) and (Hi). It is clear that Condition {ii){b) of Lemma 7 with g 7^ or 

Condition of Lemma 7 with to 7^ entail that T is conjugated with 

a piecewise affine map. Then Condition {ii){a) of Lemma 7 (with g 7^ 0) or 
Condition (i)(a) of Lemma 7 (with to / 0) cannot hold for a QMQ system 
satisfying UNI (a). 

Finally, assume UNI (a) and Condition {i){a) of Lemma 7 with to = 0. Then 
relations of Condition (i)(6) of Lemma 7 taken at the fixed points of h 
imply that exp[iToc(/i)] = 1 for all h eTi. This is only possible if c is lattice 
of span L and tq is a multiple of 27r/L. □ 
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3.6. Final bounds for the Dirichlet series S{s, w). With Relation 1)2. 17() 
which relates the Dirichlet series S{2s, w) to the quasi-inverse of the operator 
^, we obtain now the expected properties for the Dirichlet series S{s, w). 
The first result is relative to the case when w is near 0, and will be useful 
in Section 4, while the second result is relative to the case when w belongs 
to a compact vertical segment, and will be useful in Section 5. 

Lemma 8. Consider one of the three algorithms G,IC,0, and a cost c of 
moderate growth. For any with < ^ < 1/5, there is ao > and, for aU 
So, with < So < ao, there are a (complex) neighborhood W of and a 
constant M" such that for any w G W', the following holds: 
(i) ^a{w) > 1 — (ao — So). 

(a) The (meromorphic) functions s i— > S{2s,w), s S{2s,w) have a single 
pole at s = a{w) in the strip \^s — 1| < ao, and this pole is simple. 

(m) max^|5(2s,u;)|,|5(2s,u;)|^ < M" max(l, , Vs, = l±ao. 

Proof. Let W be a complex neighborhood of in which a{w) from Propo- 
sition 0(7) is well-defined. Each vertical line 3fts = o" is split into three: 

JVear the real axis. For {s,w) in a (complex) neighborhood A of (1,0), 
Proposition 0(5) gives a decomposition Jls,w = A(s, w) Ps,w + ^s,w where 
Ps,w is rank-one and the spectral radius of N^^^ is < 9. It is easy to see 
that the (1, l)-norm of (I — N^^io)"^ is bounded by some M2 on A. Since 
(t(0) = 1, taking a smaller neighborhood W2 C W of 0, there are a2 > 
and t2 > so that the following set is contained in A 

A2 := {{s,w);w e >V2,|5Rs- 1| < a2, \Qs\ < t2} ■ 

For (sjw) G A2, the quasi-inverse of H^^^ satisfies 

(3.35) (/ - H,,^)-i = ^ ^^'f^ . P.,^ + (/ - N,,^)-i . 

1 — A(s, w) 

It has as only singularities in A2 a simple pole at each point (s = a{'w),'w), 
with residue the nonzero operator 

(3-36) n{w) := \ P^(^),^ . 

Next, note that since A'^(1,0) / we have i := inf^^ | | > 0. Fix a 

with < a < a2. Up to taking a smaller W2, we have \^a{w) — 1| < ^a/2 
for w £ W2- Thus for = 1 ± a, \t\ < t2, and w G W2, the dominant 
eigenvalue satisfies 

\X{s,w)-l\ >£\s-a{w)\ >£\^s-^a{w)\ >y, 

and, by we have ||(I - Us^n])~^\\i,i < M2/a. 

Compact region. Suppose that t 7^ 0. We first prove that d{l, SpHi+j^^o) > 
0. Fixing p > p, the spectrum SpHi+j^^o decomposes into two parts 
= SpHi+i(_o n {|A] < p} and 5+ = SpHi+jj^o n {|A] > p}. Proposi- 
tion 0(1) implies that is a finite set of eigenvalues of finite multiplicity 
and Proposition 1 proves that 1 does not belong to S^. Then ^(1,5^"^) > 0. 
On the other hand, d{l, 5j~) > 1 — p, and finally d{l, Sp Hi+j^^o) > 0. Then, 
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by perturbation theory of finite parts of the spectrum, there exist 03 > 0, 
/? > and a complex neighborhood W3 of such that the distance between 
1 and the spectrum of H^^^, is at least P on the compact set 

^3 ■.= {{s,w); w£W3,\^s-l\ <a3, ^2 < \t\ < to} , 

where Iq = 1/ from Theorem 2. Thus {s,w) 1— > (/ — Hs^t^,)"^ is analytic 
on the compact set ^3 and its (1, l)-norm is bounded by some M3 there. 
Domain \'^s\ > l/p'^. Consider Si x Wi from Theorem 2, for our fixed ^ > 0. 
There exist ai g]0,(t(0) — fJo] and a complex neighborhood Wi oi w = Q 
such that any pair in {(s, u)); w G Wi, — 1| < qi} has its real part (cr, u) 
in Si X 1^1. 

Choose first ao '■= min(a2i 03, next W4 so that '^a{w) > 1 — (ao — So) 
on W4, and finally W := nl=iWj- Taking M' := max(M2/ao, M3, Mi), 
with Ml as in Theorem 2, we obtain the claim for S{s,'w). For S{s,w), we 
apply H2.18() and use that ({s) is bounded on vertical strips near 5Rs = 2. □ 

Lemma 9. Consider one of the three algorithms Q,IC,0, and a lattice cost 
c of moderate growth with span L. Consider some ^ with < ^ < 1/5, 
and V > 0. Then, there are 71 > 0, and a constant Q3 such that, for any 
r G M, |t| e [v,7r/L], 

(i) s I— > S{2s,iT) and s 1— > S{2s,iT) are analytic in the strip |5Rs — 1| < 71. 

(ii) max^|S'(2s,ir)|, \S{2s,iT)\^ < Q3max(l, |t|«) ,Vs,5Rs = 1 ±71. 

Proof. Fix < ^ < 1/5. By UNI, Theorem 2 gives a > and Qi so that 
for 1 9s I > l//0^, |Ks — 1| < a, and arbitrary real r, 

(3.37) \\{I--lis,ir)-%,t<Qi-\'^s\^, 

Suppose that \t\ < l/p^ and |r| G [u,7r/L]. We first prove, as in Lemma 
8, that d(l, Sp Hi+j( jT-) > 0. Proposition 0(1) implies that the spectrum 
SpHi+jf^jT- decomposes into two parts S^^ = SpHi+j^^jT- n {|A] < p} and 
5^"*"^ = SpHi+jf^j^ 1^ {l-^] > pji where S^^ is a finite set of eigenvalues 
of finite multiplicity. Proposition 1 gives that 1 does not belong to 5^"*"^. 
Then d{\,S^^) > 0. On the other hand, d{l,Sf^) > I — p, and finally 
ti(l, Sp Hi+it^iT-) > 0. Then, by perturbation theory of finite parts of the 
spectrum, there are < 71 < a and /? > such that the distance be- 
tween 1 and the spectrum of Ho-+it,jr is at least f3 on the compact set 
k — 1| < 71, 1^1 < l/p^, \t\ G [v,tt/L]. Then s 1— > (/ — 'H.s,iT)~^ is an- 
alytic there, and there is Q2 so that 

(3.38) ||(I-Hi±^,+it,i,)-i||i,i <g2,V|r| G [t;,7r/L], V|t| < l/p2. □ 

4. Limit Gaussian distribution for costs of moderate growth. 

In this section we prove our Central Limit Theorem, Theorem 3. We first 
explain our use of Perron's formula (§4.1). We next introduce in §4.2 a 
smoothed model endowed with probability Ptv- For this smoothed model, 
Lemma 10 allows us to deduce from bounds on the bounds on which 
entail quasi-power estimates. Theorem gives asymptotic normality for the 
smoothed model, with asymptotic estimates for its expectation and variance 
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(§4.3, in particular Lemma 12). A comparison of the uniform and smoothed 
distributions of (Lemma 14 in §4.4) finahy yields Theorem 3. 

4.1. Using Perron's formula. Choose < < 1/5. Lemma 8 specifies 
ao £]0, 1/2] and W'. We fix w G W', the Dirichlet series S{s,w) being 
viewed as functions of s. Consider the strip S{w) limited by the two vertical 
lines 3?s = 1 — ao and 5Rs = 1 + ao- By Lemma 8, this strip contains 
s = a{w) G M as unique (simple) pole of S(2s,w). In the rectangle U{w) 
defined by the strip S^w) and the two horizontal lines SJs = ±U, Cauchy's 
residue theorem provides 

where E{w) = F^.^^) o R(i(;)[l](0), with Il{w) the residue operator from 
(|3.36j) . Note in particular that 

(4.2) E{0) = -^FiPi[l](0) = -^Fi[/i](0) / 0. 

We now let U tend to oo. By Lemma 8, the integral on the leftmost vertical 
line ?fis = 1 — ao exists and satisfies 

(with a 0-term which is uniform for w G W and T — > oo), while the integrals 
on the horizontal lines ofU{w) tend to zero for U oo. Finally, Perron's 



formula (|2.2fl|) with D = 1 -\- ao gives the contribution from the rightmost 
vertical side, so that 

(4.3) xi>^^{T) = —-^^^ r^2a{w)+i h+o/'r-^So' 

^ ' ^ ' a{w){2a{w) + 1) I V 

with E{0) 7^ and a uniform 0-term with respect to w G W , as T — > cxd. 

4.2. Smoothed costs and transfer of estimates. To exploit the esti- 
mates (|4.3|) on the Cesaro sums ^^{N), we introduce an auxiliary model, 
the smoothed model. 

Associate to some nonnegative function T ^{T), with e(T) < 1, the 
probabilistic models (r2Ar(e), PAr(e)) as follows: For any integer N, set 
07v(e) = ^N] next, choose uniformly an integer Q between A^— [A^e(A^)J and 
N, and draw uniformly an element {u,v) of Oq. Slightly abusing language, 
we refer to the function C in the model (17 at (e), Pat (e)) as the "smoothed 
cost." The cumulative value of exp[?i;(7] for Pjv(e)) is 



N 

(4.4) ^-(^y-=jN7m ^ ^''^^"^ 



=N-lNe{N)\ n<Q 

SO that the moment generating function of the smoothed cost is just 
(4.5) EA.[exp(u;C)] = |4^. 
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Note that can be expressed as a function of via 

(4.6) ^UN) = Y]^^^ [*-(^) - ^-(^ - ■ 

Now, to transfer the bound (|4.H|) for into a bound for we shah appeal 
to a result that is often used in number theory contexts. 

Lemma 10. Let W he a complex neighborhood of 0, and let Cn{w) he 
a sequence of complex-valued functions on W. Assume that ^^{T) := 
Yln<T ^ri{w){T - n) satisfies 

^UT) = F^{T) [1 + (G(r))] , T ^ oo , 

with a 0-error term which is uniform for w G W, where Fyj{T) = B{w)T"'^'^^ 
and B{w), a{w) are hounded holomorphic functions such that '^a{w) > 1, 
B{w) ^ on W. Assume further that G{T) tends to for T ^ oo and is of 
moderate variation, i.e., there exists K so that \G{cT)\ < K\G{T)\ for any 
c with 1/2 < c < 2. 

Then, if G{T)-^ = 0{T) for T ^ oo, we have 
(4.7) 

yTQ^Ty,^^ - - LTG(r)i/2j)] = F'^{T) [l + 0{G{Tf') 

where the 0-term is uniform with respect to w £ W. 
Proof. We first show (without using G(T)"i = 0(r)) that for T ^ oo 
(4.8) 



*^(r)-^'^(T-rG(r)i/2) =F^{T)G{T)y^ l + OiG{T) 



1/2 



Nl/2 



1 

T 

(without the integer parts) with a uniform 0-term for w £ W. Consider 
some sequence e(r) which tends to 0. The estimate of ^'^^.(r) and the 
assumption on G entail 

^ -[^'^(T)-^^(r-re(r))] 



Te{T) 



^ {F^{T) - F^{T - Te{T))) + {F^{T)G{T)) 



Te{T) ' ' "'^ ' Te{T) 

tr^ . - ^r^^Km 1 F.UT)GiT)\ 

^w[^ ) 



^ ^ ^ V'^^^ F^TY Te{T) FUT) J_ 
Then our assumptions on Fyj{T^ and. (i{w^ imply 

Fi{T) = e {T-'F^T)) , F^(r) = e iT-'F^{T)) , 

with a uniform [recall A{T) = Q{B{T)) as T — > oo means that there is 
an absolute constant G > so that A{T) < CB{T) and A{T) > CB{T) as 
T ^ oo]. Therefore, we obtain (|4.8|) by taking 

Te{T) := ( ^'"^^}^^^ \'^ = 9 {TG{Tf'^ 



To finish, remark that the difference between (|4.8|) divided by G{T)^^'^ and 
(|4.7j) can be split into two terms of order 

FUT) • o (^^^) = KiT) • o (G(r)V2) _ □ 
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We now apply Lemma 10 to the smoothed costs. From 1)4. 3p . (|4.6() . upon 
setting G{T) = T-^^o (recall < So < 1/2), we find 



^(N) = :^MAr2^(«') \i + o( N^^A 
a{w) L \ / 



where £'(0) ^ and the 0-term is uniform (as tends to oo) when w varies 

in a sufficiently small neighborhood of 0. For w = 0, one has 

(4.9) 

^o{N) = E{0)N^\i + o(n~^A] , with ^(0) = -"^^"^^^^^^°^ 



A'(l) 

Finally, by (|4.5j) . we obtain: 

Lemma 11. [Quasi-powers for smoothed cost] Let < 70 < ao with oq 
from Lemma 8. The moment generating function of the smoothed cost 
corresponding to e{N) = N~"'° satisEes 

(4.10) EN[exp{wC)] = — ^i!^i^iv2(-{-)-{0)) [1 + (Ar-7o)] ^ 

with E{w) from i)4. J|) and a uniform O-term when N ^ 00 and w is near 0. 

4.3. Asymptotic Gaussian law for the smoothed cost. In view of 
Lemma 11 which provides a quasi-power expression for the moment gener- 
ating function of the smoothed cost, we may apply Theorem to Cn = C\^^ 
and JIat = Qat, Pat = Pat with Pn ■= logiV, kn = A^"'^°, and 

E{w) 



(4.11) U{w) = 2{a{w) - a{0)), V{w) = log 



E{0)a{w) 



We see from the above that function a{w) (which solves X{a{w),'w) = 1) 
plays a central role. The next lemma expresses in particular that U"{0) = 
2(t"(0) > 0, so that Theorem can be applied: 

Lemma 12. [Computation of constants.] In the QA4Q setting we put 
U{w) = 2{a{w) — o-(O)). Then, for the constant cost c = 1 (recaUing A{s) = 
log A(,s) from Proposition 0), one has 

2 2 9^"(^) 

(4.12) ,,= t;'(0) = ^ = ^, ,^:=U"m = ^>«. 

More generally, recalling Jl{c) and5^{c) from Theorem 1, and setting x{c) = 
A'^'.„^(l, 0), we have 

Mc) := C/'(0) = /X • /2(c) , 

d\c) := U"i0) = -il\c)-6^ + fi-P{c)+fi^flic)-xic)>0. 

In particular, U{w) (or equivalently a{w) ) is absolutely convex at 0. 

Proof. Let us begin by the case when c = 1. Then a is defined by A(a{w)) = 
—w. Therefore 

1 „, . a'{w)A"{a{w)) 

a [w) = - , a [w) = — — , 

A'{a{w)) A'^{a{w)) 
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and, recalling that A'(l) < (from Proposition 0) and A"(l) > (from 
Proposition 1), we get 

,/ N 2 , 2A"(1) 

Let us now study the general case. Taking the derivatives of the relation 
A{a{w),w) = 0, one obtains 

= a'{w)A'g{a{w),w) + A'^{a{w),w) , 

= a"{w)K{a{w),w)+a'\w)A%{a{w),w) 

+2a'iw)A'Ua{w),w)+Al,{aiw),w) . 

Remark that {a{0),0) = (1,0) and the derivatives with respect to s satisfy 

A'M^),w)\^=o = A'(l), A'^,{a{w),w)\^=o = A"(l) . 

Thus, setting L{w) := A(l + a'{0)w,w), we find ct"(0) = p^L"(0). 
Since Proposition 1 implies that L"(0) > 0, we get the strict positivity of 
(52(c), as claimed. Finally, C/'(0) and U"{0) are 

(4.13) /.(c) = 2a'(0) = , 

(A^A^ ^ 9 'Vm 2Ag(l,0)A"(l) , 4A;(1,0)A;-^(1,0) , 2A:;^,(1,0) 

(4.14) 5 (c) = 2a (0) = |^,(i)|3 + + |A'(1)| ' 

and, using 1)4.12(1 as well as Theorem 1, they may be expressed as functions 
of /i, 5, /i(c), 5(c), and x(c)- Note that x(c) = <5(c) = for constant c. □ 

By Lemma 12, Theorem applies and provides the following result: 

Lemma 13. Let < 70 < cto, with cto the constant in Lemma 8. The 
smoothed cost C associated to e{N) = N~'^° has an asymptoticaUy Gaussian 
distribution, with speed of convergence 0{l/\/\ogN). Moreover, 

(4.15) ljv[C] = C/'(0)logA^ + y'(0) + 0(iV-T«), 



(4.16) YN[C] = U"{0)logN + V"{0) + O{N~"">), 

with V defined from the residue function E{w) in i|4.1)) through If4.11\) . and 
U{w), U'iO), and U"{0) > as in Lemma 12. 

Moreover, for each fixed k > 3, there is a polynomial of degree exactly k 
with coefficients depending on the derivatives of order at most k at of U 
and V, so that the moment of order k satisGes 

(4.17) Er.iC"] = P,{logN) + O , 



with a O-term uniform in k. 
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4.4. Proof of Theorem 3. The next lemma will allow us to deduce The- 
orem 3 from Lemmata 12 and 13: 

Lemma 14. Suppose that limAr^oo e(iV) = with e(iV)"^ = 0{N/logN). 
Then the distance between the distributions Pat and PAr(e) on ri^v is 
0{e{N)). The same holds for the distance between Ptv and its smoothed 
version. 



Proof. First recall that there is ii' > so that 



for all {u,v) G Oat. The estimate provided in §1.3.4 gives K = S/ir'^ for 
g and O, and K = 3/{4tt^) for /C. Also, |f^Ar| = KN^{1 + 0{N-^)). 
Denote by A^' := N — [Ne{N)\ . Then Jljv decomposes into the ordinary 
subset On := Q^' and the exceptional subset Sjy '■= ^n\ On, with 

^n{£n) = 0{ ]=0{e{N)). 



Now, the probability ¥n{u,v) satisfies, for any {u,v) € On, 

N , N 

Em' 



¥n{u,v) = P,(n,7;)-P[g = (7] = Y 

q=N' 

Then, using N = 0{N'), the relations 

N N 
q=N' ' 



1 



+ We{N)\ 



q=N' 



g=N' 



1 + 



log A 
A' 



1 



A' 



A^ + 1 



(4.18) 



1 + 



log A 
A 



(1+LA6(A)J)|^ 



give, for any {u,v) £ On, 

- A 

Fn{u,v) = ¥n{u,v)— (1 + 



log A 
A 



FNiu,v)il + OieiN))). 



Finally, for any A C ^In, the difference | ¥n{A) — ¥n{A)\ is less than 
\FN{AnON)-^N{AnON)\ + \'PN{An£N)-^N{An£N)\ =0(e(A)). □ 

Proof of Theorem 3. For e(A) = N~'^°, Lemmata 13 and 14 imply that 
the asymptotic distribution of C is Gaussian, with a speed of convergence 
in 0(l/-v/log A). This proves Theorem 3(a) for any 7 < 70 < ao- 

To get claim (5) of Theorem 3, we use the fact that for each of our three 
algorithms, the maximal number of steps P{u,v) on ^n is bounded by 
Kq log A. This is a classical result due to the fact that the three Euclidean 
divisions we consider are of the form v = mu + er with r < v/2. For 
{u,v) £ flN, all the quotients rrii are at most A, and the moderate growth 
assumption then implies that 



sup{C''{u,v);{u,v) G Un} < K{logN) 



2k 
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for some K and all N, and finally, Lemma 14 gives 



E 



n[C'']-En[C'']\ = O ((logiV)2*^7V-To) . 

Applying Lemma 12, we get the expression for EjvfC^] and for the variance 
in Theorem 3(6), ending the proof of Theorem 3. □ 

5. Local Limit Theorem for lattice costs. 

This section is devoted to proving Theorem 4. Consider an algorithm with 

associated dynamics in the good class and satisfying UNI. Let c be a lattice 
cost, of span L > 0, and of moderate growth. By Lemma 14, it suffices to 
consider the smoothed model Fn associated to e(A^) = N~^, for 7 satisfying 
< 7 < ao (ao given by Lemma 8) which is to be fixed later. By Lemma 
11, we have a quasi-power expression for Kn[cxp{wC)] for small \w\. 

We first consider costs c with span L = 1, so that 

Our starting point is the relation 

¥j,[C{u, v)=e] = ^J^ e--^ • EAr[e-^] dr . 

Since we are looking for an LLT result, the convenient scale of the problem 
is n := log N, and we set 

^ = Px{n) := lJi{c)n + 5{c)x^fn, q^in) := [j)x{n)] . 

(Here [■] denotes the nearest integer function.) We consider 

In := 27rVlogAr • P^[C(u, v) = g^(log TV)] 

/ + TT 
exp[-irg^(ra)] • Ejv[e*^'^] dr . 
-TT 

Our strategy is to decompose the integration interval [— tt, +ir] into a neigh- 
borhood of zero, [— ^",1;], and its complement |r| G (f, tt]. This gives rise to 
two integrals In^ and In^ . 

5.1. Estimate of the integral /^'^^ We begin with and a neighborhood 
W' as provided by Lemma 8. We assume that [—v, +v] is contained in W' 
and study 

/(°) ■.= V^ [ \M-irqx{n)] ■EN[e'^'^]dT . 

J—v 

Lemma 11 grants us a quasi-power expression for EAr[exp(iTC)]. Recalling 
that ij,{c) = 2a' (0), we see that the integrand can be written as a product 
Qnfni where the first factor is 

^"^^^ E{Q)alT) ■ + ■ ^^P[^^(^^-(") - ; 

and the second factor is 

fniT) '■= exp[2n((T(ir) — 1 — iTa'{G))\ ■ exp[— zra;(5(c)\/n] . 
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The function z ^ (t(z) — 1 — z(t'(0) has a saddle point at z = 0, and the 
main useful properties of gnij)-,fnij) are 

(5.2) \gnij) - 1| « |r| + exp[-n7], \fn{j)\ < exp[-nr^(5o] , 
for |r| < V, where v is taken sufficiently small, so that 

(5.3) 6o :=inf{|5Rc7"(ir)|;|r| < t;} > . 



logn 
6on 



1/2 



We first prove that the dominant part of ^ is 

exp[-iTqx{n)] ■EN[e''^^]dT, with r„ := 
Properties ()5.2() entail that 

sup{|5n(-r)/„(r)|; t„ < |r| < z;} « exp[-nr^5o] « ^ , 

and this implies that \In^ — Jn^\ is 0{l/^/n). 
We now define the following quantities, 

/•Tn 

jj^' := \/n / exp[— ira;5(c)\/n] • exp[2n(cj(ir) — 1 — ircj'(O))] dr . 



r-Tny/n 



exp[-iT95(c)2; - -^'^ 5^ {c)] , 



:= r e^^[-m{c)x-UH\c)]d{} = ^e-^"'\ 
J-00 2 6{c) 

and will show that they satisfy \Jn^ — Jn~^^^\ = 0{l/^/n) for 1 < i < 3. Let 
us check this line by line. 

From Jn^ to Jn \ Properties 1)5.2(1 imply that 



- 4°) I « r |r| expi-nr^^oldr + e-"^ T" exp[ 



-nr^(5o]dT , 



and the two terms are 0(l/n). 

From Jn^ to Jn\ There is first a change of variable ?? = T^/n. Remark that 

(log n)^/^ \ 



tends to for n ^ 00. Furthermore, since d {c) = 2a" {0), 
2n 



n \ n 



2 V 

so that the difference between the two integrals satisfies 

|JW - « ^ ^3exp[-^2^]rf^ = o(^)- 

From Jn^ to Jn^. The difference between the two integrals satisfies 

UP - 4^)| < / eM-^'^-^W = o(^). 
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(3) 

The last equality for JA is obtained by "completing" the square. 
Finally, the first part In'' of the integral /„ satisfies 

(5.4) /(o)=^£_^+o(^ 

5.2. Estimate of the integral when r is not close to 0. We now con- 
sider the integral (|5.1() outside of the neighborhood of zero. 

Lemma 15. Assume that c is a lattice cost of moderate growth with span 
L. For each v > 0, there are Q < oo and 70 > so that, for the smoothed 
cost corresponding to e{N~'^) with small enough 7, we have 

|l7v[exp(irC7)]| < QN'^" ,V|r| G [v,7:/L] . 



Proof. We apply Lemma 9 with Perron's formula, keeping the notations 
of Lemma 9. Consider, as in Section 4.1, Perron's formula (|2.2()|) with D = 

1+71, 

'^^r{T)= [ S{2s,iT)^-^ds 
J5Rs=l+7i S{2s + 1) 

In the rectangle U defined by the strip 3f?s = 1 it 71 and the horizontal 
lines = itC/, the function s 1— > 5(25, it) is analytic (from Lemma 9) and 
Cauchy's residue theorem gives 

2s+l 



/ Si2s,w)—, -ds = Q. 

Ju ^ 's{2s + l) 



Furthermore, using once more Lemma 9, 

/ Si2s,w)- -ds = 0(r3-27i) . 

We let U tend to 00. The integrals on the horizontal lines ollA tend to zero, 
Perron's formula gives the right side, and finally 

(5.5) ^ir{T) = 0{T^-'^^^). 

Now, if 7 < 271, setting 70 = 271—7, the smoothed quantity (^ir^N) satisfies 

(5.6) = 0(iV^-i)0(iV3-2Ti) = O(iV2~70), 
Furthermore, with (jOJ, ^o(iV) = 0(iV^) . Therefore 

E7v[exp(irC)] = = 0{N-^°) . □ 



We deduce from Lemma 15 that In^ is 0(e 



-7170^ 
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5.3. End of proof of Theorem 4. Finally, we take oq and a neighborhood 
W' as granted by Lemma 8. We assume that [—v, +v] is contained in W' 
and is sufficiently small to entail validity of ()5.3() . Then, v is fixed and we 
choose 7 so that 7 is less than 271 (from Lemma 9) and less than oq (from 
Lemma 8). 

Now, with (|5.4|) and Lemma 11, Theorem 4 is proven in the case of span 
L = 1. Remark that if c has span L, then d := c/L has span 1, with 
/i(c) = Lfj,{d) and 5{c) = L6{d). This proves Theorem 4 in the general 
case. □ 

5.4. The nonlattice case. We say that a nonlattice digit-cost is shifted 
lattice if there exist L G and a number Lq, with Lq/L irrational, so that 
(c — Lq)/L is integer- valued. The number Lq is called the shift. The span 
of c is the largest possible L. Then, the digit cost c can be written as the 
sum of two costs c = Lq + d, where d is lattice with span L, and the total 
cost C equals c = LqP + D. It is then possible to apply the previous result 
to D. Furthermore, since the two A functions relative to c and d are related 
by Ac(s, w) = wLq + Ad{s, w), one has 



fi{d) = ^(c) - fiLo, 5\d) = 6\c) + Ll6^ - Lo[x(c)/i2 + 2^5^] ^ 



and, finally: 

For the algorithms Q, fC, O and for any shifted lattice cost c of span L and 
shift Lq, which is of moderate growth, one has, for any x G M 



with a 0-term uniform in x. 

The case when c ^ is neither lattice nor shifted lattice is essentially dif- 
ferent, the noncompact situation requiring Dolgopyat-type estimates in w. 
The generic nonlattice situation will be considered in a forthcoming paper. 



This article has presented a unified approach to a large body of results 
relative to Euclid's algorithm and its major (fast) variants. We recapitulate 
here some of its consequences, relations to the vast existing literature on 
the subject, as well as extensions and open questions. It should be stressed 
that most of the improvements can eventually be traced to the existence of 
pole-free strips for Dirichlet series of cost parameters, a fact that precisely 
devolves from our extension of Dolgopyat's estimates to continued fraction 
systems. 

First our methods lead to extremely precise estimates of the moments of 
costs, and in particular the number of steps, a much studied subject. Our 
estimates, when specialized to the mean number of steps, yield 



Pjv {u,v);--< C{u, v)-Lq- P{u, v) - fi{d) log N - 6{d)xy^logN < - 




6. Conclusion. 



(6.1) 



En[C]= nlogn + 7] + 0{N-^); 
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see above Theorem 3, Parts (6) and (c). In the case of the standard algo- 
rithm, this covers the original estimates of Dixon and Heilbronn in 1969- 
1970 (the main term), as well as Porter's 1975 improvement (the second 
term, ij), while providing the right shape of the error term {0{N~^)), for 
which Porter further showed that one could take 7 = ^ — e in the case of the 
standard algorithm. We refer the reader to the accounts by Knuth [HOI and 
Finch ^Sl for more material on this classical topic. Our formula (|6.1j) also 
extends Rieger's analyses (first published around 1980, see [1211101) of the 
centered algorithm and the odd algorithm. Note that, in our perspective, 
the second-order constant rj comes out as a spectral quantity. It is an open 
problem to derive an explicit form starting from our expressions (e.g., such 
a form involving C'(2) is known under the name of Porter's constant in the 
standard case). 

In sharp contrast to the mean value case, variances do not seem to be 
amenable to elementary methods. The first-order estimate has been given 
by Hensley (1994) in the paper that first brought functional analysis to 
the field. Our formula, 

(6.2) YN[C]=6Hogn + 6i + 0{N-^), 

stated in Theorem 3 can be viewed as a precise form of Hensley's estimate 
relative to the standard algorithm, one that also extends to the odd and 
centered algorithms. (Incidentally, the quantity 5, called Hensley's con- 
stant ^HI) is not recognized to be related to classical constants of analysis, 
though it is now known to be polynomial-time computable thanks to a re- 
cent study of Lhote |S2j; the nature of 61 is even more obscure.) Note that 
the complex-analytic properties of the moment generating functions pro- 
vided by the functional-analytic methods furnish similarly strong estimates 
for moments of arbitrarily high order (our Theorem 3), a fact which also 
appears to be new. 

Regarding distributional results, several points are worthy of note. Dixon 
in his 1970 paper had already obtained exponential tail bounds and these 
were further improved, albeit in a still fairly qualitative way by Hensley in 
his 1994 study (see his Theorem 1 in j24j ) . Our approach gives much more 
precise information on the distribution, as we now explain. 
For simplicity, let us specialize once more the discussion to the number of 
steps. First, regarding the central region of the Local Limit Theorem (The- 
orem 4), the nature of the error terms obtained {0{N^^^'')) and the fact 
that the saddle point method lends itself to the derivation of full asymptotic 
expansions (see, e.g., Henrici's book [221) entail the existence of a full asymp- 
totic expansion associated with the Gaussian approximation of Theorem 4, 
namely, for a computable numeric sequence {cj}, 

(the argument of Pat being that of Theorem 4). This expansion is valid for x 
in any compact set of M. 

Regarding large deviations. Lemma 11 (quasi-powers for smoothed costs) 
implies that E,iy[ex.p{wC)] is of the form of a quasi-power, provided w 
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stays in a small enough fixed neighbourhood of 0. By a standard argu- 
ment (originally due to Cramer and adapted to the quasi-powers framework 
by Hwang |251 1261 07] ^ and Lemma 14, this implies the existence of a large 
deviation rate function, that is, of a function I{y) such that, for the right 
tail corresponding to y > 0, there holds 

(6.4) hm -^logP^[C > (l + y)^logiV] = -I{y), 

N^oo log iV 

with I{y) defined on an interval [0, 5] for some 6 > 0. In other words, 
the probability of exceeding the mean by y// log is exponentially small in 
the scale of the problem, being of the form exp[— /(y) log A^] = N^^^y\ A 
simplified version of ()6.4() is then: the probability of observing a value at z 
or more standard deviations from the mean is bounded by a quantity that 
decays like e~'-''^^ (Hensley's notations j24j ) . Analogous properties naturally 
hold for the left tail. Cramer's technique of "shifting the mean" and use of 
the resulting "shifted" quasi-powers lead in fact to very precise quantitative 
versions of 1)6. 4|) in the style of Hwang — thereby providing optimal forms of 
tail inequalities d la Dixon-Hensley. 

Similar properties hold for other costs measures, like the ones detailed in the 
introduction. Of particular interest is the statistics of the number of digits 
assuming a particular value, for which estimates parallel to 1)6. If) . (|6.2() . (|6.3|) . 
and (|6.4|) are also seen to hold. For instance, the frequency of occurrence of 
digit m in the expansion of a rational number has mean 

~ log2 ( 1 + , ^ I OA I ' 

and it exhibits Gaussian fluctuations. It is again unknown to us, which of 
the involved constant (beyond the mean value factor) may be related to 
classical constants of analysis. The spectral forms of our Lemma 12 may at 
least provide a starting point for such investigations. 

A major challenge is to derive distributional information that are so to speak 
"superlocal" . By this, we mean the problem of estimating the behaviour of 
the number of steps and other cost measures over fractions of denominator 
exactly equal to A^, i.e., rationals belonging to \ J^at-i. In view of 
what is known in the average-case |18|l3f)j. we expect arithmetical properties 
of N to come into the picture. The analytical difficulty in this case stems 
from the fact that a further level of "unsmoothing" (how to go from to 
Oat \ r^Tv-i?) would be required. 

Further works. Our dynamical approach provides a general framework, 
where it seems possible to answer other questions about distributional anal- 
ysis. For instance, all the questions that we solve for rational trajectories 
can be asked for the periodic trajectories. A periodic trajectory is produced 
by a quadratic number, and the reference parameter is related to the length 
of the geodesies on the modular surface. In a forthcoming paper, and fol- 
lowing the same principles as in j371 1381 HT] , we prove that the distribution 
of costs on periodic trajectories can be studied in a similar way as here, 
replacing the quasi-inverse (/ — H^^tu)"^ by det(/ — H^^^). 
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Open problems. We also ask various questions about Euclidean algo- 
rithms: for instance, what happens for other Euclidean algorithms of the 
Fast Class (in particular for the Binary algorithm [Hl|^)? The extension of 
our results to cost functions that are still "small" but take into account the 
boolean cost (also known as the bit complexity) of each arithmetic operation 
is on our agenda. Note that an average-case analysis is already known to 
be possible via operator techniques |21lini- On another register, the exten- 
sion to "large" costs is likely to lead us to the realm of stable laws: see for 
instance Gouezel and Vardi's works \21\ I5()j for occurrences of these laws in 
continued fraction related matters. 
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